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Application of Satellite Data in Observational and Theoretical 
Studies of the Evolving Structure of Baroclinic Waves 


1 . Introduction 

During the contract period we completed a variety of 
observational and theoretical studies aimed at elucidating the 
relationship between satellite measurements of cloud and 
radiation and the evolution of transient and stationary features 
of the circulation in middle latitudes. 

In broad outline, our efforts can be divided into four main 
categories: (1) the application of satellite radiation data in 

global weather analysis, (2) studies of blocking weather 
patterns, (3) studies of finite amplitude baroclinic waves 
including hydrologic processes, and (4) the development of a low- 
order dynamical systems model for the nonlinear evolution of 
stationary and transient planetary waves. A listing of the 
papers published with the support of this grant, and the 
complementary USRA grant NAS8-34010, follows. 

I. APPLICATION OF SATELLITE RADIATION DATA IN GLOBAL WEATHER 

ANALYSIS 

1) Hansen, A.R. and R.L. Nagle, 1984: Estimates of the 

generation of available potential energy by infrared 
radiation. Mon. Weath. Rev. , 112, 1370-1377. 

2) Nagle, R.L., 1984: Climatological aspects of the 

outgoing longwave radiation field of the northern 
hemisphere mid- latitudes. M.S. Thesis , Yale 
University, pp. 59. 

II. STUDIES OF BLOCKING WEATHER PATTERNS 

3) Hansen, A.R. and A. Sutera, 1984: A comparison of the 

spectral energy and enstrophy budgets of blocking 
versus nonblocking periods. Tellus , 36A, 53-63. 
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4) Hansen, A.R. and A. Sutera, 1984: Observational 

aspects of the predictability of atmospheric blocking. 
Predictability of Fluid Motions , (G. Holloway and B.J. 
West, eds.) American Inst. Physics, 441-447. 

5) Benzi, R., A.R. Hansen and A. Sutera, 1984: On 

stochastic perturbation of simple blocking models. 
Quart. J. Roy. Met. Soc. . 110, 393-409. 

6) Brown, P.S., A.R. Hansen, and J.P. Pandolfo, 1986: 
Circulation regime-dependent nonlinear interactions 
during northern hemisphere winter. J. Atmos. Sci., 43, 
476-485. 

III. STUDIES OF FINITE AMPLITUDE BAROCLINIC WAVES INCLUDING 

HYDROLOGIC PROCESSES 

7) Tang, C.-M. and G.H. Fichtl, 1983: The role of latent 

heat release in baroclinic waves - without b-effect. 

J. Atmos . Sci. , 40, 53-72. 

8) Tang, C.-M. and G.H. Fichtl, 1984: Non-quasi- 

geostrophic effects in baroclinic waves with latent 
heat release. Atmos . Sci. , 41, 1498-1512. 

9) Saltzman, B. and C.-M. Tang, 1985: The effect of 

finite-amplitude baroclinic waves on passive, low- 
level, atmospheric constituents, with applications to 
comma cloud evolution. Tellus , 37A, 41-55. 

10) Tang, C.-M. and B. Saltzman, 1985: Comma cloud 

development in a two-layer baroclinic model with 
vertical water vapor flux at mid- level. Papers in 
Meteorol . Res. , 8, 45-52. 

2. Summary of Research Results 

In two major applications of an 8-year record of satellite- 
derived outgoing longwave radiation (OLR) we have (1) computed 
the generation of available potential energy due to infrared 
radiative fluxes and (2) examined the extent to which these data 
can provide signatures of high and low frequency weather 
phenomena including blocking. Using conventional temperature 
data along with the satellite OLR data it has been estimated by 
Hansen and Nagle (1984) that the global pattern of infrared 
radiation causes a destruction of both zonal and eddy available 
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potential energy of magnitudes G z ~ - 5.0W m“ 2 and G £ ~ - 0.6 w 
m“ 2 respectively. The destruction in the eddies is almost 
entirely due to the very long stationary disturbances of 
wavenumbers one and two. In R.L. Nagle's Masters thesis it is 
demonstrated that the OLR data contains the signatures of 
intermittent blocking phenomena embedded in higher frequency 
baroclinic wave activity that can be correlated with observed 
longwave amplified states. 

Significant contributions to forming the observational 
foundation for a theory of blocking are provided in the sequence 
of papers by A.R. Hansen and collaborators. For example, Hansen 
and Sutera (1984a, b) show the differences in the energy and 
enstrophy cascades in blocking and non-blocking situations, 
establishing that pronounced upscale flow of both of these 
quantities, from intermediate to planetary scales, occurs during 
blocking episodes. The upscale flux of enstrophy, in particular, 
suggests that the persistence of blocking periods may be due to 
reduced dissipation of the large scale circulation and therefore 
entail some above normal predictability. The observational 
results (Brown, Hansen and Pandolfo 1986) also indicate that 
smaller scale transient eddies play an important, regime- 
dependent, role in interactions with atmospheric circulations on 
the scale of blocking. 

In a theoretical study (Benzi, Hansen, and Sutera, 1984) 
some of the effects of these smaller scale, higher-frequency, 
eddies are modelled as both a stochastic white and red noise 
forcing of a bistable orographical ly induced stationary flow of 
the type studied by Charney and Devore. An important new time 
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constant thus emerges, namely the 'exit time' for reversal 
between the two attractor domains of the steady flow representing 
blocking and normal states. These exit times were found to range 
widely as a function of the zonal forcing, and in the red noise 
case it is found that the bimodal nature of the barotropic model 
can be obscured depending on the prescribed variance of the 
perturbations and the model parameters. 

Our third area of work concerned the development of a 
minimal model of finite amplitude baroclinic waves and of their 
satellite signatures. These signatures usually take the form of 
evolving cloud patterns that, in the mature stage, often have the 
well-known "comma" structure. In Saltzman and Tang (1985a) we 
discuss the dynamical processes involved in such comma cloud 
formation, as well as the more general question of the 
redistribution of a low-level passive tracer in a growing 
baroclinic wave. An analytical solution has been obtained 
demonstrating the complex nongeostrophic flow pattern involved in 
(i) the redistribution of constituents trapped in the lower 
layers of a finite amplitude baroclinic wave, and (ii) in the 
formation of the typical humidity and cloud distributions in such 
a wave. This model is extended in Tang and Saltzman (1985b) to 
include the vertical flux of moisture across the middle (e.g., 

500 mb) level to and from the upper level. The results show that 
there is a tendency for the upper-level maximum mixing ratio to 
shift to the east relative to the lower-level maximum mixing 
ratio. This is consistent with the eastward tilt with height of 
the common form cloud in the "cold-conveyor" cross-section and 
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the poleward tilt of the cloud in the "warm-conveyor" cross- 
section as observed by Carlson (1980). 

3. Recommendations and Conclusions 

The emerging results from our contract studies, and 
complementary work done by other scientists associated with our 
efforts (e.g., Sutera, Speranza, Benzi, Hansen at Control Data 
Corporation) , provide a new outlook on the problem of low- 
frequency weather variability, one that offers some hope that 
increased predictability may be in the offing if we can properly 
treat the longwave bimodality that has been established. We have 
been pursuing this new opening vigorously under our new contract 
now in effect (NAS8-36356). In particular we have been 
developing a generalized low-order model capable of yielding 
solutions exhibiting bimodality in planetary wave structure along 
with a unimodal probability distribution in the zonal wind field. 
In its ultimate form this model will include explicit reference 
to the hydrologic processes that determine the cloud and 
radiative fields measureable from satellite platforms. 

At the same time we have been pursuing observational studies 
of time-series satellite data using mathematically-sophisticated 
methods of attractor analysis to determine the minimum level of 
complexity needed to model the data and to determine the extent 
to which rapid transitions associated with bimodality can be 
identified in the data. 

Results from these new approaches are not reported here, but 
are discussed in the Quarterly Reports of contract NAS8-36356, and 
in a new series of papers published, in press, or in preparation. 
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APPENDIX 

Reprints of papers published with support of NASA contract 
NAS8-34903 . 
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OBSERVATIONAL ASPECTS OF THE PREDICTABILITY 
OF ATMOSPHERIC BLOCKING 

Anthony R. Hansen 

Yale University, New Haven, CT 06511 
Alfonso Sutera 

Center for the Environment and Man, Inc., Hartford, CT 06120 

ABSTRACT 

Enstrophy and kinetic energy flux functions calculated from 
observed data from the winters of 1976*77 and 1978*79 are examined 
to compare the nonlinear cascading properties of blocking and non* 
blocking periods. During 1978*79 a more pronounced upscale cascade 
of kinetic energy as veil as an upscale cascade of enstrophy from 
intermediate to planetary-scale wavenumbers was found during block 
periods, as compared to nonblocking. During 1976*77, similar re- 
sults appeared for a case of Rex blocking but not for the persis- 
tent, greatly amplified planetary wave pattern in January and 
February 1977 (as Identified by Charney, et al. 1 ). 

The predictability time 2 based on the enstrophy flux function 
showed an Increased predictability for the January- February 1977 
event but no significant differences between the other blocking 
cases compared to the nonblocking sample. However, the reversal of 
the low wavenumber enstrophy cascade during blocking does suggest 
that blocking may be more persistent (due to reduced dissipation of 
the large-scale circulation) and therefore more predictable. Some 
possible implications for theoretical modelling are discussed. 

INTRODUCTION 

In the present study, we will attempt to Identify any system- 
atic differences in the energy and enstrophy cascading statistics 
of blocking events compared to nonblocking periods using data from 
two recent winters (1976*77 and 1978-79). It is generally accepted 
that the existence of persistent features such as blocking highs 
will allow better medium range forecast skill in numerical weather 
prediction models. Recently, Bengtsson 3 has shown from numerical 
simulations that the atmosphere is In general more predictable dur- 
ing blocking events than during nonblocking periods. Ve are inter- 
ested to see if any diagnostic evidence exists for the greater pre- 
dictability of blocking found by Bengtsson by qualitatively esti- 
mating a predictability time based on two dimensional turbulence 
theory. Although our sample Is fairly small, some interesting 
features appear. 

This report is part of a more complete set of diagnostic 
calculations to appear in Tellus 4 . 

0094 - 243X/84/1060441 - 07 $3.00 Copyright 1984 American Institute of Physics 
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DATA AND PROCEDURES 

The data used in this study are the horizontal wind (u, v) from 
the twice daily operational analysis of the National Meteorological 
Center for the winters of 1976-77 and 1978-79. Fourier coefficients 
of u and v were computed for every 2.5* of latitude at the 10 man- 
datory levels in the troposphere with the wavenumber expansion trun- 
cated after wavenumber 18. 

Following Steinberg et al. 5 , the kinetic energy and enstrophy 
flux functions can be defined as 

6F (m) 

IS — " - Sc (“I 11 * 1 ) (i) 

and 

C £ («n|n.l) (2) 

assuming F R (0) - F £ (0) - 0. (m|n,l) and Cj (m|n,l) represent 

the gain in wavenumber m kinetic energy and enstrophy, respectively, 
due to nonlinear, triad interactions with all possible combinations 
of wavenumbers n and 1. These equations were integrated from 30*N 
to 80*N and from 1000 mb to 100 rib. 

Our primary goal is to study persistent, large-scale blocking 
events compared to predominantly zonal circulations which may in- 
clude short duration, low amplitude, small-scale features. There- 
fore, we restrict ourselves to cases of stationary or slowly propa- 
gating ridges where the departure of the 500 mb height from the 
zonal mean, averaged over every 2.5° of latitude from 55°-80*N, 
exceeded 250 m for 7 days or more. Any observation not falling in a 
blocking period was included in the nonblocking sample. 

In addition, a case of very large negative height departures 
occurred during January and February 1977. Ve will use the dates 
given by Gharney, et al. 1 for 2 periods of persistent height depar- 
tures from the climatological mean in our diagnostic calculation. 
Because of the existence of large amplitude features of one type or 
another throughout the 1976-77 winter, no nonblocking days from this 
winter are included in our nonblocking sample. A summary of the 
blocking and nonblocking days for the 2 winters is given in Table L 

Comparison of the periods of large, persistent height depar- 
tures with synoptic charts Indicates that these periods satisfy (for 
the most part) the conventional, subjective definition of blocking. 

RESULTS 

The results for the 1978-79 winter are given in Figure 1. Nor- 
mally, enstrophy is cascaded from the longest to the shortest wave- 
lengths 5 . This characteristic also appears in our nonblocking 
sample (Fig. la, dashed line). However, the absence of an enstrophy 
flux out of the lowest 5 wavenumbers during blocking and the change 
in sign of the enstrophy flux function near wavenumber 2 is evident 
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Table I Tabulation of blocking and nonblocking days from the 
1978-79 winter and the 1976-77 winter determined 
from persistent height departures 

1978-79 


blocking days (90 observations) 


0000 

GMT 

i 

Dec. 

- 

1200 

GMT 

7 

Dec. 

0000 

GMT 

20 

Dec. 

- 

1200 

GMT 

27 

Dec. 

1200 

GMT 

29 

Dec. 

- 

0000 

GMT 

10 

Jan. 

0000 

GMT 

14 

Jan. 

- 

1200 

GMT 

26 

Jan. 

0000 

GMT 

16 

Feb. 

- 

0000 

GMT 

22 

Feb. 

nonblocking 

days 

(85 observations) 

0000 

GMT 

8 

Dec. 

• 

1200 

GMT 

19 

Dec. 

0000 

GMT 

28 

Dec. 

- 

1200 

GMT 

28 

Dec. 

1200 

GMT 

10 

Jan. 

- 

1200 

GMT 

13 

Jan. 

0000 

GMT 

27 

Jan. 

- 

1200 

GMT 

15 Feb. 

1200 

GMT 

22 

Feb. 

- 

1200 

GMT 

28 

Feb. 


1976-77 

blocking days (36 observations) 

1200 GMT 11 Dec. - 0000 GMT 29 Dec. 

Chamey et al.* negative anomalies (76 observations) 

30 December - 17 January 
2 February - 22 February 


in Fig. la (solid line). Also, a much larger flux of kinetic energy 
from the Intermediate to the long waves occurred during blocking 
as can be seen clearly in Fig. lb. 

A comparison of the energy and ens trophy flux functions for 
the December 1976 case of blocking and the negative anomalies of 
January and February 1977 are given in Fig. 2. The December 1976 
event exhibited a very striking reversal of the low wavenumber en- 
strophy flux function (Fig. 2a), and the same enhanced upscale 
kinetic energy flux as the 1978-79 blocking cases (Fig. 2b). The 
similarity in the behavior of Rex-type blocking in the 2 winters 
lends support to the significance of this result. 

The January-February 1977 ens trophy flux (Fig. 2a, dashed line) 
is more like the 1978-79 nonblocking sample. However, notice the 
much lower ens trophy flux function for the January-February 1977 
case in a nearly constant range from roughly wavenumbers 8 to 15. 

The enstrophy flux function in the same range for the December 1976 
blocking case is virtually identical to the 1978-79 blocking and 
nonblocking samples. The kinetic energy flux for the 1977 case is 
dominated by exchanges amongst the lowest wavenumbers (Fig. 2b, 
dashed line) . 

It is well known that the rate of error growth in numerical 
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F E 

4 . '78-'79 



1 2 3 4 5 B 7 0 10 15 10 

wav«numb«r 


. 8 r 



1 2 3 4 5 6 7 0 10 15 ie 

wav«numb«r 


Fig. 1. The nonlinear flux functions for a) enstrophy, 
F £ and b) kinetic energy, F R for blocking (solid line) 

and nonblocking (dashed line) for the 1978-79 winter. 
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Fig. 2. The nonlinear flux functions for a) ens trophy 
and b) kinetic energy for the December 1976 blocking 
case (solid line) and the January-February 1977 
negative anomalies (dashed line). 
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weather prediction models due to small-scale Inaccuracies in the 
initialization can be related to the rate of downscale enstrophy 
cascade in a constant enstrophy flux inertial subrange of a two- 
dimensional turbulent fluid 6 . The characteristic time of this pro- 
cess can be called the predictability time and is inversely propor- 
tional to the cube root of the constant rate of enstrophy transfer 
to higher wavenumbers in a constant enstrophy flux inertial sub- 
range 2 * 6 . The fact that the atmosphere at the larger scales behaves 
predominantly like a two-dimensional fluid in which the kinetic 
energy spectrum obeys a -3 power law 7 * 8 * 9 suggests that a qualita- 
tive estimate of the predictability time for blocking compared to 
nonblocking situations based on the enstrophy flux function may be 
made. The estimate can only be done in a qualitative way because 
the numerical value of the enstrophy flux function, F^, in the con- 
stant range is dependent upon the truncation used 6 . The reason for 
this sensitivity is that finite amount of enstrophy is lost when the 
wavenumber expansion of the wind field is truncated at wavenumber 18. 
This will cause aliasing problems in the nonlinear enstrophy* inter- 
action in our truncated dataset that will cause changes in the shape 
of the enstrophy flux function at high wavenumbers. Unfortunately, 
the accuracy of our data at wavenumbers higher than 18 is poor so 
this problem cannot be alleviated. However, we will assume the values 
of the flux function that we have are qualitatively correct for com- 
parison purposes. 

As noted earlier, the enstrophy flux function in the constant 
range is nearly identical in both blocking and nonblocking situations 
[Figs, la and 2a (solid line)] suggesting no greater predictability 
for Rex- type blocking. However, the January-February 1977 low 
anomaly event exhibited a roughly 1/3 smaller enstrophy flux rate 
(Fig. 2a, dashed line). This result allows us to speculate that the 
persistent low anomalies in 1977 were inherently more predictable 
than either the nonblocking or Rex blocking cases. 

An explanation of the greater predictability of the Rex-type 
blocking found by Bengtsson* might be found in the reversal of low 
wavenumber enstrophy flux function, and the absence of an enstrophy 
flux out of the lowest 5 wavenumbers during blocking as opposed to 
the normal case. This may indicate that the destruction of large- 
scale vorticlty by down-scale cascade and eventual dissipation is 
greatly reduced or eliminated during blocking events. As a result, 
blocking patterns may be more persistent and therefore more predict- 
able. In addition the enstronhy cascade reversal may be an indica- 
tion of intermittent behavior* 0 in the large-scale flow requiring 
no particular low wavenumber Instability. 


CONCLUDING REMARKS 


Because of the evidently strong interaction between barodlnlc 
cyclone-scale waves and planetary-scale waves during blocking 4 * 11 , it 
would appear that theoretical models should incorporate the effects of 
transient, barodlnlc eddies. An effort in this direction has recent- 
ly been put forward by Benzl et al 12 . Using a low-order, barotropic 
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model 13 , they suggest that the inclusion of baroclinic activity as a 
parameterized forcing of the mean flow (as a red process), acts to 
mask the role of any underlying orographic instability that otherwise 
might be present. 
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Summary 

In this paper we discuss some problems connected with the stochastic dynamics generated by the effect 
of small random perturbations acting on prototype equations governing large-scale flows. We focus on the 
features of the process which occur when the main instability is the result of orographic forcing on simple 
barotropic motion. In particular, we study a three-mode truncation of barotropic flow over topography when 
a stochastic white forcing acts on the system. We found analytical estimates of the exit times from the attraction 
domains of the steady solutions of the model. The exit times were found to range widely as a function of the 
zonal forcing. 

In addition, we included a representation of the large-scale forcing in terms of a red process. We found 
that as a function of the variance of the perturbing process and the model's parameters, the bimodal nature 
of the barotropic model can be obscured. 


1. Introduction 

A most striking feature of complex phenomena such as atmospheric motion is the 
enormous number of degrees of freedom cooperating to determine a particular observed 
pattern. On the other hand, if we are concerned with long enough time scales, we may 
observe that only a few of the many degrees of freedom play a relevant role in determining 
the long-term patterns. It is, therefore, relevant to design a model where the large-scale, 
long-term behaviour is predicted by governing equations applied to a finite number of 
components of the system and to parametrize the short-term, smaller scales as a small, 
random perturbing process. 

Schematically, the governing equation for the n-active components assumes the 
form 

i*b(x) + W (1) 

where x is an n-dimensional vector, b is an n -component vector field, Wisa n-dimensional 
stochastic forcing, and the over-dot denotes a time derivative. 

It is customary to measure the strength of W by its variance £(W 2 ) = e(e might be 
time- and x -dependent), where £ is the mathematical expectation operator. We recall 
that the total effect of W must be small, therefore e will be assumed to be a small 
parameter; i.e. small in relation to certain properties associated with the nature of b. 

If we consider W as a white process, (1) is a system of stochastic differential 
equations. Since Einstein’s application of them to the theory of Brownian motion, such 
systems have provided a good framework for describing deterministic motion in a random 
environment, in particular the macroscopic variability. 

Inspection of (1) reveals that the stochastic perturbations force the deterministic 
(e* 0) motion, changing instantaneously the local (in time) behaviour of the system. 
However, the most important effect of their action is that, in the long term, they may 
dramatically change the averaged motion. These facts were discussed by Sutera (1980) 
“Present address: European Centre for Medium Range Weather Forecasts, Reading. 


393 



394 


R. BENZ1, A. R. HANSEN and A. SUTERA 


in the framework of the Lorenz (1963) model (which remains a remarkable test-bed 
model for investigating features that may be important in the atmosphere). They were 
further considered by Moritz and Sutera (1981). Here we wish to stress a few quantitative 
facts that were neglected in the previously mentioned papers (which addressed different 
aspects of the same problem, such as the role of metastable chaos in prototype models). 

If e, as we said, is the variance of the noise, the large departures from the purely 
deterministic motion occur only on the time scale 0(e _1 ) = T t . Hence, within the time 
interval (0, t < T,) no appreciable differences between deterministic and stochastic 
behaviour would be expected. What will happen for times longer than T t depends again 
on b, i.e. on the laws governing the evolution of the averaged quantities (that is, the 
average taken over an ensemble of realizations of W). Since departures from this average 
depend on the work done by the random forcing against the deterministic forcing, it is 
of predominant interest to estimate this quantity. This paper is, therefore, dedicated to 
the calculation of the work done by stochastic perturbation of simple low-order models 
of orographically forced barotropic flows. We shall see that T t , for these models, ranges 
widely as a function of the zonal forcing and that the equations obtained by considering 
steady forcing are insufficient for the purpose of a ‘realistic’ description of the phenomena 
studied. 

We suggest inclusion of a time variability of the forcing by considering it to be a red 
process with a given amplitude and decorrelation time. We find that, in this case, the 
probabilistic nature of the barotropic model may be so much changed as to obscure its 
bimodal structure. 

2. The models 

In this paper we confine our study to the behaviouT of a simple barotropic oro- 
graphically forced flow in a /5 channel and exclude models on the sphere such as those 
of Kallen (1981, 1982) and Legras and Ghil (1983). Moreover, we will concentrate on 
highly truncated variants of these models, excluding studies where more modes are 
considered, by comparing the behaviour of two different models: those of Chamey- 
Devore (1979) (hereafter CDV), and Egger (1982) (hereafter E). 

We first write the governing equation. Let H be the mean height of a homogeneous 
fluid on a plane. Let y and x be the northward and eastward directions, respectively. 
The fluid is confined in a /3 channel by zonal walls at y - 0 and y-nL. The lower 
boundary elevation is described by z - h(x,y). Let r; be a perturbation of the upper free 
surface. If the motion is quasi-geostrophic then, following Chamey (1974), the governing 
equation is given by conservation of potential vorticity 

j t F 2 V- V'/A J ) + /{V'.(?V+ (V2tf)/o +&>)}= -KVXv-V) (2) 

where V 2 is the Laplacian operator, y - grj// 0 , g is the gravitational acceleration gravity, 
fo = 2£2sin <po, = 45 e, , A 2 = gH/fl, p = (2Qcos<po)/a where a is the earth’s radius 
and Q the angular velocity of the earth’s rotation, and J(a,b) = (da/dx)(db/dy) - 
(da/dy)(db/dx). 

The right-hand side of (2) contains the effect of Ekman dissipation (~KV 2 y) and 
a fictitious vorticity source (KV 2 y*), simulating a large-scale forcing presumably of 
baroclinic nature. In reviewing the nature of low-order models, Lorenz (1982) described 
the procedure to be adopted for obtaining such models by starting from the continuum 
equations describing the fields. It consists of considering a set of orthonormal functions 
satisfying the same boundary conditions that the fields satisfy with the additional con- 
straint that if /„ and f m are any two functions of this set then (/*/«, ) * where the 
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brackets denote a horizontal average. The expansion is performed by assuming 

ip(x,y,t) = '2ip m (i)f mt . (3) 

mn 

CDV and E chose f m as the set of eigenfunctions of the Laplacian operator satisfying 
the condition of no normal flow at the boundaries. With this choice, the set is readily 
constructed with trigonometric functions. Moreover, the same set is employed in expand- 
ing the stationary forcing ip* and the topography h. If ip^ and h m are the coefficients 
of the expansion in ip* and h, by simple substitution (2) is transformed into 

dip/dt = b(ip,ip*,b), (4) 

where ip, ip* and h are vectors in the infinite dimensional space spanned by f m and b 
is a vector field on it. The passage from (2) to (4) is straightforward, but the statement 
(3) is not because it implies that whatever/,* may be, is a steady solution of 

(2). To our knowledge there are few cases in which (2) has known steady solutions of 
the form (3), and they are restricted to very specific physical situations. However, it can 
be demonstrated that in (2) if K is very large, then asymptotically (in K) we may neglect 
the Jacobian term and hence the above-mentioned choice of /«* is properly made. 

Nevertheless, neglecting the previously discussed difficulty, (2) is transformed to 
the form (4). As an additional assumption to make any analysis of (4) tractable, we must 
require that all but a few components of the vectors ip, ip* and b be zero for all time. 
By employing this assumption we reduce our problem to studying a finite system of 
ordinary differential equations. In CDV and E only three orthonormal functions are 
considered, reducing the infinite system to one with only three equations. The functions 
are: 

gA * V2. cos(y/L), g* = 2 cos(nx/L) sin(y/L), g L = 2 sin(nx/L) sinCy/L). 

The topography is represented by the expansion h(x,y) - h 0i cos(nx/L)sin(y/L). 

In this paper, the stationary forcing is assumed to be ip*(x,y) = ip% cos(y/jL). 
We use dimensionless time scale »■»./, space scale x (or y) = x (or y)/L, and h = 
h/H. The expansion for the dimensionless ip is given by 

V = ipAgA + $KgK + rpLgL ; 

and the corresponding low-order model is 


<pA * -K(ipA - $A ) + *01 <f>L 

(5a) 

-b nl ipL- RipK 

(5b) 

ipL * Ript- *„i ipA 

(5c) 


where 

*oi - (WO + A- J ))(V2H) h ni - ( W(n 2 + 1 + A - 2 ))(h 0i /2H) 

b n i = rp A ~ Pm - (8 V2/3 Jt)(n 3 /(n 2 + 1)) 

Pm = ( L/a)(n/(n 2 + 1)) cot <po y*i m (8 V2/3« )n . 

Hereafter we choose n *= 2, R = 10" 2 and L/a * i, as in CDV. The form of system (5) 
is a feature common to E. 

Since we force the barotropic flow only in its zonal component, it is legitimate to 
consider stochastic perturbations only on that component. We notice that in this case 
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the stochastic forcing acts on the total angular momentum only and this assumption 
might be considered unrealistic. However, recent observations (Barnes etal. 1983) seem 
to indicate that, besides its well-known seasonal oscillations, the total angular momentum 
is randomly fluctuating with an amplitude even larger than the one here assumed. 

Of course, a more realistic assumption would consist of stochastically perturbing 
the wave equations (i.e. V'x and \p L ) also; for the sake of simplicity we restrict ourselves 
in forcing only the mean flow. By choosing a Wiener process as the perturbing random 
force, system (S) describes a diffusion process. In the next section we shall study it in 
detail for the CDV model. The results carry over identically to the E case. Henceforth, 
we will drop the overbar notation for dimensionless quantities. 

3. Diffusion against the flow 

It is known that the model (5) exhibits multiple equilibria, specifically that the 
steady-state equations possess three solutions x b , x„,, and x z , where x ■ (ipA,Vic, VO. 
for selected values of the mountain height h and driving forcing xpX ■ Moreover, x b and 
x z are asymptotically stable solutions, and x m is an unstable one. If D\ and D 2 are the 
attraction domains of xt, and x z respectively, then an initial condition chosen in D\ or 
D 2 will approach x b or x z respectively (see Fig. 1). 



Figure I. Schematic diagram of the phase space for the CDV model. The symbols are explained in the text. 

The effect of stochastic forcing on time scales T t will be to drive trajectories starting 
from any initial conditions to points outside the domain of attraction to which those 
initial conditions belong. How do the trajectories behave for small e? 

For xo in Du x(r) will roll down into a neighbourhood of x b along the deterministic 
trajectory. The trajectories then stay in a small neighbourhood of Xt,, making excursions 
to other points of D\ from time to time and returning. Finally, after a long time, T„ the 
trajectory deviates so far from xt, that it goes over into D 2 and rolls down into a 
neighbourhood of x z . Then the process will begin again. This process is completely 
described by the joint probability distribution P. Recalling that the noise acts only on 
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y> A , the Fokker-Planck equation for P is 

BP(v>a, V>k, n>L, t)/bt = ied 2 P/dy>\ - div(b. P) (6) 

where b = {-K(y A - rp A ) + AoiV^ - b nl y> L , -Ky L + b„iyt K - h„iy> A }, and div 

is the three-dimensional divergence operator defined in the (y> A , y> K , VO space. E solves 
(6) by numerical integration. 

In this paper we will try to estimate P in a different fashion. First we notice that 
two time scales are present in our system, namely l/K and 1/e. If l/K < 1/e, in the sense 
discussed in appendix A, the asymptotic behaviour of our problem is approximated by 
the single stochastic differential equation 

dy> A = ~( dV/dy> A )dt + e^dW (7) 

where V(y A ) is given by 
- 

V(vu) = - J F(y> A )dy> A 

c 

= “ J [{-(V^ - Va)K 2 - blj(ip A - IPX) - h Hl h 0i y> A }/(2K + bh/K )] . dy> A . (8) 

In appendix A we propose a way to obtain (7) by employing a technique resembling the 
diffusion limit of the theory of stochastic differential equations (Schuss 1980, pp. 132- 
134). We are aware that to solve our problem, other asymptotic techniques are available 
(Schuss 1980), and we have employed them finding the same result. Since our purpose 
is to use our result for a wide range of the parameters (leaving, however, K and e 
unchanged), we have tested our approximation numerically and the results are presented 
in appendix B. It appears from these results that the approximation is valid. The situation 
described by the approximation (7) is identical to the one discussed by Sutera (1981) and 
the steady joint probability distribution is given by 

P(Va) = N' exp(-2V(v>„)/f) (9) 

where N' is a normalization constant. 

However, for (7), in the limit e—* 0, a very important statistical quantity can be 
calculated. Let us consider the process r b (x b ,x) defined by 

thfo, x) = {if x(0) is in £>i , x(r b ) is not in Z>, }, (10a) 

lb being the exit time of the process starting in D x . Analogously, t*(x„ x) is defined by 

r,(x 2 , x) = {if x(0) is in D 2 , x(r,) is not in D 2 }. (10b) 

As discussed in Sutera (1981) for the case of three steady states of which two are stable 
.'solutions and the third is an unstable one as in (7), the following estimations of the 
expectation value of the exit times hold: 

£{r b (v>$,vu)} * £{T b (^, • V-(V'S)l“*exp(2A 1 V/£) (lla) 

EM V'i, V > A )} « EMv'a, Vtf)) - *\V”(v* A ) . V'(y>X) |-» exp(2A 2 V/e) (lib) 

where y A , y A , ip A are the blocking, zonal, and intermediate equilibria, respectively, V" 
denotes the second derivative of V with respect to xp A calculated at the arguments, 
A\V- IV(V'a) -V(v»”)| and A 2 V= |V(v>Ji) -V(yj)| . (Equations (11a, b) hold for 
the limit e—* 6. Benard LeGras (private communication) has calculated the next order 
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(in 1/e) correction for larger e. We have not employed his correction in this paper in 
order to limit its length. We thank LeGras for his communication.) 

Equation (9), together with (11a, b), allow exploration of the parameter space on 
which V(y A ) depends. Here the main differences between this work and E are found. 
In concluding this section we are amused to remark that the main approximation 
employed (i.e. Eq. (7)) has been assumed valid in the same spirit in which the orthonormal 
basis has been chosen. In fact, the equation 

/( v>, V v + fth/H + fy) « - JC(V V - V V) ( 12 ) 

in the limit 1/K — * 0, converges to 

V 2 (y- V>*) = 0. (13) 

Therefore, eigenfunctions of the Laplacian are also steady solutions of (2) in this limit. 


4. Exploring the parameter space 

In this section we employ the analytical formulas developed in the previous section 
(and verified in appendix B) to study the parameter space on which the models of CDV 
and E depend. We choose to study the behaviour of the exit time and probability 
distributions as functions of h 0 i/H and the zonal forcing (yZoruo for E case). Hence 
all the other parameters are left unchanged from the ones adopted by CDV and E. 


TABLE 1. Barrier heights and exit times for the Charney-DeVore model as 

A FUNCTION OF THE ZONAL FORCING, Wa 

(a) hoi/H - 0*05 


Wa 

(n.d.) 

V(l)-V(2) 

(xl0- s ) 

V(2)-V(3) 

(xlO- 5 ) 

r (block) 
(days) 

r (zonal) 
(days) 

■p 

2*17 

0*004 

76*6 

180*1 


1-76 

0-61 

18-1 

43*6 


1-48 

1*79 

10*8 

49*8 

0*16 

1-24 

3*43 

7*6 

79*2 

018 

105 

5*51 

5*8 

159*5 

0*20 

0*88 

8*02 

4*6 

395*6 

0*22 

0*73 

11-00 

3*9 

1*2 x 10 3 

0*24 

0*61 

14*30 

3*4 

4*3 x HP 

0*26 

0*49 

18*10 

3*0 

1*8 x HP 

0*28 

0*40 

22*20 

2*8 

9*6 x HP 

0*30 

0*31 

26*80 

2*7 

5*8 x 10 s 


(b) hm/H - 0*10 


0*14 

11*10 

0*13 

1-3 x 10 s 

81-9 


10*40 

1*15 

532*7 

59*6 


9*83 



84*7 

0*20 

9*30 

4*91 

189*3 

164*0 

0*22 

8*83 

7*52 

128*7 

402*9 

0*24 

8*40 

10*60 

92*0 

1-2 x 10 3 

0 26 

8*01 

14*10 

68*1 

4-4 x 10 3 

0-28 

7*64 

18*00 

51*8 

1-9 x 10* 

0*30 

7*30 

22*40 

40*4 

1-1 x 10 3 

0-32 

6-98 

27*20 

32-1 

6-7 x 10 3 


The indices 1, 2, 3 refer to blocking, orographic instability, and zonal steady states. 
t (block) and r (zonal) are the times for the blocking-to-zonal and zonal-to-blocking 
transitions, respectively. 
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TABLE 2. Barrier heights and exit times for the Egger model as a function of 

THE ZONAL FORCING. U Q 

(a) hoi/H = 0 05 


uZ 

V(l)- V(2) 

V(2) - V(3) 

r (block) 

r (zonal) 

(ms' 1 ) 

(xl0- 5 m : s-') 

(xlO-’mV 1 ) 

(days) 

(days) 

7*0 

0*94 

0*13 

18*6 

57*1 

80 

0*80 

0*47 

12*1 

43*2 

90 

0*69 . 

0*95 

9*0 

41*8 

10*0 

0*59 

1*55 

7*4 

45*9 

11*0 

0*52 

2*27 

6*2 

55*1 

12*0 

0-44 

3*10 

5*5 

70*7 

13*0 

0*38 

4*04 

5*0 

96*4 

14*0 

0*32 

5*09 

4*6 

139*3 

15-0 

0*27 

6*25 

4*3 

212*5 

(b) hoi/H - 

010 




10*0 

6*01 

0*02 

216*8 

159*7 

11*0 

5*67 

0*32 

92*7 

65*6 

12*0 

5*39 

0*84 

62*2 

58-1 

13*0 

5*15 

1 52 

46 1 

62*4 

14-0 

4*92 

2*34 

35*6 

74*1 

15*0 

4*72 

3*29 

28*9 

96*4 

16*0 

4*53 

4*36 

23*6 

133*8 

17*0 

4*36 

5*50 

20*0 

196*0 

18-0 

4*19 

6*64 

17*0 

289*8 


The indices 1* 2, 3 refer to blocking, orographic instability, and zonal steady states. 
t (block) and r (zonal) are the times for the blocking-to-zona) and zonal-to-blocking 
transitions, respectively. 






Figure 2. The CDV model potential, V(v^). for hoi/H m 0 0$ as a function of the zonal forcing yZ * 
(a) yZ * 012; (b) yZ * 0*14; (c) yZ * 0*18; (d) yZ m 0*20, (e) yZ * 0*24. Psia denotes y A . 
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Figure 3. The CDV model potential. V(vu), for WW = 0-10 as a function of the zonal forcing yX 
(a) v; * 014; (b) yX - 0-18; (c) yX - 0-20; (d) VM* - 0-24; (e) vtf - 0-26. 



Figure 4. The CDV model potential V(y A ), for hoi/W-0-20 as a function of the zonal forcing yX . 

(a) yX - 0-22; (b) yX - 0-30; (c) yX - 0-40. 


Moreover, the noise level in dimensionless units has been fixed following E as t - 
5 x 10 -5 . The results are reported in Tables 1(a) and (b) for the CDV case (where a 
transformation into natural units has been performed for the exit times), and Tables 
2(a) and (b) for the E case. Figures 2 to 4 contain the potential V(ip A ) for the CDV 
cases (the E case does not present remarkable differences). It can be noticed that the • 
basin of attraction of the zonal state is a very shallow one. The feature has been captured 
by employing the term {2K + (b nl b nl )/K} as a weighting factor. 

Recently, Dole and Gordon (1983) analysed the statistical nature of anomalies in 
the atmosphere and found that potentials of the nature presented here would fit the 
data. However, according to Dole and Gordon a shallow basin would describe a ‘blocked’ 
circulation rather than the zonal one we found. 
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5. The inclusion of red noise 

A first look at the tables presented in the previous section shows that although there 
are combinations of parameters which make the models quite attractive, slight changes 
in the zonal forcing are able to change interesting quantities such as exit times in a 
dramatic fashion. 

Barotropic models have been employed traditionally to gain a first insight into the 
phenomena under investigation. The need for a better modelling of the zonal forcing 
is a peculiar characteristic of these models and since the fictitious source of vorticity 
corresponds to no particular physical process in the governing equation , we may speculate 
about its nature. 

A possible assumption may be formulated by considering this forcing to be a 
.stochastic process obeying a stochastic equation. In this section, we shall consider the 
case when the zonal forcing, y> A , is itself the solution of a stochastic process. We will 
model Va by employing a forcing with a red spectrum. This case has relevance in our 
problem because in the atmosphere this forcing may be associated with barodinic 
processes. There is some observational evidence that blocking events are initiated and 
influenced by a forcing from transient, barodinic, cyclone-scale waves (see Hansen and 
Chen 1982; Hansen and Sutera 1984). 

For simplidty we postulate that 

= + (e k /K)dW (14) 

with f as previously defined and fa a constant. Let us assume that y obeys 

dy * - Amydt + (e 1 Ao>) , dW. (15) 

Hence the process y will be an Omstein-Uhlenbeck process (see, for example, Sutera 
1981) with decorrelation time 1/Aa>. 

We shall study the following problem: 

d\f>A = — {dV'fvv* ; v A ,y)/dy A )dt + e*dW (16a) 

dy * -A wydt + (tiA ot^dW. (16b) 

If Aco ► e,, then on time scales Ofef 1 ) the process y may be considered 6-correIated; 
i.e. 

E\y] * 0 (17a) 

E[y(>)ym - (e,/2Ao>)6(/), (17b) 

and it is Gaussian. Hence, (16a, b) reduce to 

dif> A ~Pdt + g(y/ A )ydt + e*dW. (18) 

The functions P and g are defined as follows: 

/= (~K\va - fa) - bli(it>A - fa) - h 0 ih ni y> A )(2K + b’,/*)' 1 (19a) 

g( Va) - (K 2 + bh)(2K - bii/K)~\ (19b) 

In (18) the stochastic forring enters in a twofold way. We have additive (f*dW) and 
multiplicative (g(V' A )y) sources. In the limit described by (17a, b), the appropriate 
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Fokker-Planck equation corresponding to (18) is given by (see Schuss 1980, p. 95) 


±r.-JL(f + 


£\ d 


at 


dy A \ 2Ao)dy 


W) 


p*l 1 


2 2Acody>A 


(R+g 2 (H>A))P 


where R - 2Awe/e t . A stationary solution of (20) is 

1 4 Aa> f 


p ’ ( exp — J rT1^7) 


( 20 ) 


( 21 ) 


We plot the probability distribution, P{rp A ), obtained from (9) where only additive 
white noise is present (e = 5 x 10' 5 ) compared with that obtained from (21) in Figs. 
5(a) and 5(b), respectively. In Fig. 5(b) we have taken = 2 x 10"\ Aoj = 0-02 (cor- 
responding to a decorrelation time of ~5 days) and £ = 5 x 10' 5 . The inclusion of both 
multiplicative and additive noise (Fig. 5(b)) has the effect of leaving the maximum' 



Figure 5. Normalized probability distribution for the parameters « 014, hm/H « 0 05 and K * 0 01. (a) 
P(Va) from (9) with e« 5 x 10~ s ; (b) P,( ) from (21) with t*5x 10' 5 , t, - 2 x 10 _< and Aw - 002; (c) 

P,( y > A ) with Ci - 2 x 10' 4 , Aw* 0 02 but with e « 0 (i.e. no additive noise); and (d) same parameters as (b) 
but where g(vvt) has its minimum on the zonal side of the instability. (This change was obtained by setting 

<K,i * 1 and is intended as an illustration only.) 
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likelihood of the blocked state unchanged while reducing the maximum likelihood of the 
zonal circulation compared with a system in which only additive white noise is present 
(Fig. 5(a)). This effect is associated with the function b„ j. In fact, the function g(y> A ) 
has a minimum where b„i - 0, namely, where y> A = {P„i/a ni ). This value of y>A lies 
between the blocking equilibrium and the orographic instability. Therefore, if only red 
noise is present (Fig. 5(c)), the process will tend to remain on the blocked side of the 
double potential well. All this is expressed by (21 ) and confirmed by numerical calculation 
(not shown). If the point where g( rp A ) is minimum fell on the zonal side of the orographic 
instability, the system would tend to remain in a zonal configuration (Fig. 5(d)) in the 
absence of additive white noise (e = 0). Extended numerical integrations of both the 
approximation (16a, b) and the full CDV mode (5) in which Va is modelled as in (14) 
yield histograms of the residence density of VU that are qualitatively similar to the 
.analytical results presented in Figs. 5(a) to (d). 

These results lead to some speculation about the nature of an atmospheric signal. 
According to the simple barotropic model proposed by CDV, we should expect that 
' such a signal would be bimodal. However, if we record this signal in a barodinic 
environment as in the atmosphere, it is not obvious that such bimodality should be 
observable unless we were able to filter the red component of the signal. Hansen and 
Sutera (1983) compared statistics for blocking and non-blocking periods and, although 
no evidence of bimodality was present in their statistics, active barodinic processes were 
present both during blocking and non-blocking periods. Thus, in the context of the 
results presented here, this barodinic activity might mask the effect of an underlying 
barotropic mechanism which would otherwise generate a bimodal signal. 

6. Concluding remarks 

We have presented a set of results concerning the use of stochastic perturbations 
in simple models of barotropic flow over large-scale topography. We found that stochastic 
forcing may drive the motion against the deterministically stable flow. We were able to 
calculate the rate at which different equilibria would be visited as a function of the 
models’ parameters. If we interpret the multiple equilibria of these models to represent 
different equilibrium states of the atmospheric drculation (e.g. low and high index 
equilibria), the following conclusions may be drawn from our calculation: 

First, the exit time calculations show that such equilibria may be so energetically 
separated that the distribution of their occurrences may be meaninglessly uneven. 
Moreover, even when their exit times are roughly in agreement with the observed 
duration of these patterns, a small change in the external, large-scale flow forcing has 
a dramatic effect. When the dynamics of the fordng are simulated by some red process, 
we noticed that the bimodality of the probability distribution may be obscured. 
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Appendix a 

In this appendix we describe the procedure used in deriving Eq. (7). The model is 


ij>A m -K(V a~ Va) + hoiy> L + e*W (A.la) 

Vfc = -K\p K - bm^ L (A. lb) 

i>L * -KVl + b„ iVk - h ml rj> A . (A.lc) 

We introduce the transformation 

V m Va) + hoiV't- (A.2) 

In these new coordinates, the system (A,la, b, c) transforms into 

VU-rj+cW (A. 3a) 


V>* = -Ky L - (Whoi){»? + K(Va ~ V'.S)} (A.3b) 

>)= —2Kt) — K\y A - y>Z)+ b nl h Ql y K - h 0 ih nl y A - KeW (A.3c) 

We notice that in (A. 3a, b, c) two time scales are involved, namely K and t. Moreover, 
K/e is a large parameter (0(10*)) using the values considered in the text. 

Setting t = a, the previous system is transformed to 

Va * (l/e)r?+ W (A.4a) 

\i> K = -(K/e)y> K - (b„i/eh 0l ){T] + K(\f> A - y>i)} (A.4b) 

17 = -(2K/e)r\ - (l/e){-K\y A - n) + b Hl h 0 m- hoi h n m)- KW (A.4c) 

where the fact W(t) ~ W(r/e) * £*W(r) has been used. We may formally integrate 
(A. 4a) to obtain (all integrals run from 0 to r) 

V'^V'o+O /e) J yds* J Wds. (A.5) 

Let 

n ■ * + y (A. 6 ) 

where 

dy/dr = - (2 K/e)y + K t) (A.7a) 

and 

y(0) * »7(0), (A.7b) 

then 

dtj/dr = dy/dr + dx/dr = -(2 K/e)y + KW + dx/dr. (A. 8 ) 

It follows that 

dx/dr = (2 K/e)y - KW - (2K/e)r 1 - 

- {l/e){-K\\t> A -it>Z) + b ni hoirt>K-h 0 ik„iV A } + KW (A.9) 
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or 

(dx/dr) = -(2 K/e)x - (l/e){ -K\y A - yX) + b n fi 01 rt, K - h 0 fi ni rp A l (A.10) 
Equation (A. 10) has the solution 

x m J exp{-(2K/e)(r- s)} x 

x [-(l/*){ -K\v>a + bnihoiVic - hoih^^ds. (A.ll) 

On the other hand, the solution of (A.7a, b) is 

y » exp{- (2K/ e)T}y(0) + \e{2K/e) txp{-(2K/e)x) J exp {(2K/e)r}W(r). (A.12) 

. Now let us recall that W is not a differential function; hence the last integral of (A.12) 
must be evaluated by parts. It follows that 

* ie(2K/e) txp{-(2K/ e)x) f exp{(2K/e)r}W(r) 

- ic(exp{- (2K/e)r) | (2 K/e) exp{(2 K/e)r}W(r) | j - 

- exp{- (2K/ e)r}(4/k 2 /c 2 ) J exp{(2 K/e)r)W(r)dr] (A.13) 

- H(2K/e)W(x) - txp{-(2K/e)z)(AK 2 le 2 ) j exp{(2K/e)r}W(r)dr]. 

Substituting (A.13), (A.12) and (A.ll) in (A. 6), and taking into account (A.7b), we 
have 

TJ = exp{-(2AT/f)T}rj(0) + (1/2 K) J (2 K/e) exp{-(2K/e)(r - s)} x 

x {K 2 (V>a - V>X) + b H ih 0i y>K - h 0 ih n} yi A }ds + (A.14) 

+ W(2K/e)W(x) - txp[-(2K/e)x)W 2 /e 2 ) f exp{(2 K/e)r}W(r)dr}. 

On the other hand, from (A.4b) it follows: 

Vk - exp{- (K/e)r}y>K(0) - 

- (1/K) J (K/e) exp{-(K/e)(T-s))(bJh 0l ){Ti+K(y>A- Va))*. (A.15) 

Now we take the limit (K/e) —* « and notice that 

lim[(AT/e) exp{-(K/c)(r - j)}] -» 6(r - s). (A.16) 

In this limit, (A.14) and (A. IS) may be written as 

tj * (l/2K){-K\y A - y>Z) + b Hi hoirj>K - (A.17a) 

V* - ” (1/ K)(b n \/ ho\){rj + K(xp A - v»5)}. (A.17b) 

* Solving for tj and substituting in (A. 5) we get 

. »,-»,«> zggjj . »» - *»'*•' » (r) , (A . I8) 

After changing back to time r, (A. 18) is equivalent to the stochastic equation used in 
the text (Eq. (7)). We remark that in (A. 18) we are left only with term O (K 2 /e), i.e. 
about two orders of magnitude smaller than K/e. 
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The crucial point of our approximation resides in assumption (A. 16). The numerical 
value of K/e is 10 3 . We believe that this is a proper asymptotic limit and the numerical 
calculation in appendix B supports our belief. The next step would be to carry *on the 
expansion of the function (K/e) txp(-(K/e)(x - s)} in powers of K/e. 

We remark that the technique here employed resembles and generalizes a way to 
demonstrate the diffusion limit for Brownian motion and also known as the Smoluchowski 
approximation (Schuss 1980 p. 132, compare the introduction of the variables x and y ). 
Many other proofs of the Smoluchowski approximation are available including one based 
upon asymptotic expansion of the Fokker-Planck equation corresponding to the process. 
We have performed such an expansion and found the same result as in (A. 18). 

Appendix b 

In this appendix we present numerical evidence about the approximation derived 
in appendix A. In Figs. 6 to 8 we present sample trajectories for the full model and for 



0 100 200 300 400 500 


Tim* (day*) 

Figure 6. Sample solution of the CDV model with - 010, hoi/H * 0 05, K - 0 01 for (a) the full mode. 

(b) the approximation. 
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T i me (days) 

Figure 7. Sample solution, as Fig. 6 except y>X “ 0-14. 


the approximation. We choose three cases. In the first (Fig. 6) the blocking solution is 
most probable. In the second (Fig. 7) the two equilibria are equally probable solutions 
and the zonal solution is most probable in the third case (Fig. 8). Statistically, we remark 
at the similarity between the two signals. For r*5x 10~ 5 , we have integrated the two 
models for a time equivalent to 5787 days counting the number of exits from one domain 
to the other one. The results are summarized in Table Bl. We also calculate the residence 
time of the process. It is defined as the fraction of time spent by the trajectory in a cell 
of the phase-space. The cross-section along the axis is presented in Fig. 9. From all 
this we may conclude that our approximation has a very good performance. 

Finally we remark that our study is independent of the use of a perturbing process 
in only one component of the system. In fact, in the case of a perturbing process projected 
on all three equations, the technique is almost the same. 



Paid 







PERTURBATION OF BLOCKING MODELS 


409 


TABLE B1 . Statistics on the exit times (in days) obtained by integrating 
THE FULL MODEL AND THE APPROXIMATION FOR 5787 DAYS 



Number of 
exits 

Average 
exit time 

Standard 

deviation 

Full model 

Zonal — ► Blocking 

130 

31-40 

43-00 

Blocking — ► Zonal 

166 

6-69 

8-00 

Approximation model 

Zonal -*■ Blocking 

122 

34-70 

45-60 

Blocking — * Zonal 

133 

10-30 

10-44 


h 0 \ /H * 0 05 and V’JS * 014 


Barnes. R. T. H., Hide, R., 
White. A. A. and 
Wilson. C. A. 

Chamey, J. G. 

1983 

1974 

Chamey, J. G. and 
DeVore, J. G. 
Dole, R. M. and 
Gordon, N. D. 

CDV 1979 
1983 

Egger, J. 

E 1982 

Hansen, A. R. and Chen, T. C. 

1982 

Hansen, A. R. and Sutera, A. 

1983 

Kallen, E, 

1981 


1982 

Legras, B. and Ghil, M. 

1983 

Lorenz, E. N. 

1963 

1982 

Moritz, R. E. and Sutera, A. 

1981 

Schuss, Z. 

1980 

Sutera, A. 

1980 


1981 
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ABSTRACT 

The redistribution of a low-level, passive constituent of the atmosphere under the influence of a 
growing barodinic wave is examined by a series of analytical calculations based on a two-level, 
highly truncated model. It is shown that a constituent confined to the lower half of the 
atmosphere, and initially homogeneous in the horizontal, will tend to achieve maximum 
concentration in the low pressure/warm sector portion of the wave and minimum concentration 
in the high pressure/cold outbreak region with sharpest gradient between the maxima and 
minima roughly coinciding with the cold front. This distribution is further accentuated if an 
initial meridional gradient of the constituent exists. 

If we assume, as a rough first approximation, that water vapor can be considered to be such a 
passive constituent, it is shown that the implied relative humidity field and cloud distribution will 
tend to evolve into the comma-type form commonly observed on satellite images of mid-latitude 
cyclone waves. Moreover, the solution is shown to replicate the complex flow regime associated 
with the comma formation, elucidating the dynamical roles of vertical motion and advection in 
the cloud evolution. 


1. Introduction 

One consequence of the motions associated with 
growing baroclinic waves in middle latitudes is the 
redistribution of any passive atmospheric con- 
stituents (or tracers) that may be present at the 
time of inception of the waves. We shall investigate 
this redistribution with the two-level, 0-channel, 
model of the evolution of finite-amplitude baro- 
clinic waves developed previously (for a review see 
Saltzman and Tang (1982). The fundamental 
question we pose and discuss is the following: 
suppose a barodinic wave forms in a field 
containing a constituent of horizontally uniform 
low-level mixing ratio e, or of zonally uniform 
meridional gradient ( dt/dy ). To what finite-am- 
plitude distribution will the constituent evolve under 
the influence of the motions generated by the 
baroclinic wave? 


Secondly, we ask whether the above inferences 
are consistent with the observed distribution of 
water vapor and cloud in a barodinic wave, the 
latter of which is often characterized by a distinct 
“comma” form. It will be assumed for this purpose 
that, as a first-order approximation, we can treat 
water vapor as such a passive constituent and that 
we can identify cloud coverage with some critical 
values of the relative humidity and vertical motion. 
Although it is clear that a full dynamical de- 
scription of these cloud and moisture phenomena 
must depend on more detailed numerical models, it 
is also of value (at least from a pedagogical 
viewpoint) to deduce the essentials of the pheno- 
mena and its dynamics with the simplest model 
capable of analytical solution. In fact, one may 
argue that only from such “minimum” physical 
models do we gain a true “understanding” of 
complex phenomena. It is toward this end that this. 
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as well as our previous studies, have been directed. 
As we shall see, using our simple model it is indeed 
possible to account qualitatively for many of the 
main features of the comma cloud dynamics as 
revealed by detailed synoptic and satellite ob- 
servations (Carlson, 1980). 


2. The model 

The symbols to be used in this study are the 
same as used previously (Saltzman and Tang, 
hereafter called S-T, 1972, 1975, 1982) and are 
defined in Appendix A. For present purposes we 
now add to this list the dependent variable £, the 
mixing ratio for an arbitrary constituent of the 
atmosphere. The full set of basic equations govern- 
ing the dependent variables y (stream function), x 
(velocity potential), co (individual pressure change), 
and care (S-T 1972), 


d dtff dco 

— K.-V(VVW-r /— + * = 0, 

ct ox op 


&a> £ d 

/—— + _ V* a> - — ( V # - V(V 2 ¥ )\ 
dp* f dp 


-fi 


C* iff 

dp dx 


V 2 


v # V 


diff 


-P= 0, 


dco 


dt dt 

— + VVc+ a >- — Y = 0, 
dt dp 

where 


( 1 ) 


( 2 ) 

(3) 

(4) 


d dco 

Af = V x * V(V* ip) + fiv + co — V 2 iff — V 2 iff 

dp dp 



Y 4 = source function for c. 


As in the previous studies we consider the 
two-layer system portrayed in Fig. 1, applied to a 
fundamental harmonic region of a £plane having a 
meridional width D and zonal extent L (re- 
presenting the wavelength of maximum baroclinic 
instability). For simplicity we assume that the 
constituent is confined only to the lower layer 
wherein its vertical mean concentration is given by 
c, and the mean production Y t is assumed to be 
zero. This would be the case, for example, if c were 
constrained to the following vertical distribution. 


c(p) 


2Pt l (p-Pi) h 
0 


( P>P2 ) 
(P<P 2 ) 


corresponding physically to the behavior of a 
constituent whose initial source is near the surface 
and whose concentration always decreases mar- 
kedly with height. The variable e, is, therefore, a 
measure of the mass of the constituent in the 
complete vertical column (or, in the special case of 
water vapor, of the "precipitable water”). 

Next, we expand each of the dependent variables 
( = (y/, x , cj, e) into a zonal mean part £ (y\ 0 and 
an eddy part £' (x, y, /), as in S-T (1975), i.e.. 




and further resolve each of these parts into an 
initial, basic, zonal field denoted by the superscript 
(0); a primary, quasigeostrophic, field denoted by 
the superscript (1) (Phillips, 1954); and a second 
order, non quasi-geostrophic, field denoted by the 


p(mb) subscript 

0 0 

250 il^il I 


500 


51 


2 


750 




3 


1000 777777777777777777777 4 

Fig . L Two-level model of baroclinic flow as repre- 
sented by the dependent variables r, *, co and e. The 
pressure levels are shown at the left and the corre- 
sponding subscript notation is shown at the right. 
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superscript (2); i.e., 
i ={<«» + {"» + 
e = 

so that, e.g., for { = c, we have, 
c = e< 0) 0’) + y 9 /) + ? w 0\ 0 + ^ 2, (*>.>\ /). 

(5) 

We can then derive systems of equations for 
each of the above components, which will differ 
from those summarized in S-T (1982) only by the 
addition of equations for c. The systems of 
equations take the following forms: 



Pi 

(III) Equations for the secondary wave fields, 
l *«>]: 




U)\" 

Pi 


B BA" 

— cj” + fij #l —— 

Bt dx 


(6c) 

B^Bif ^ 
dx By p 2 

(6d) 


(11) Equations for the forced zonal state 


sff «? 2 d><*> 

7^r 


/>2 


(7a) 



Where double signs appear in these equations, 
the upper sign refers to level I (250 mb) and the 
lower sign to level 3 (750 mb); S is assumed to be a 
horizontally uniform value of the static stability 
representative of conditions prevailing near the 
wave center (i.e., near y * JZ>, and A f B , and C 
represent non-homogeneous forcing due to the 
lower order fields, the expressions for which are 
given in S-T (1972, 1975) and are restated in 
Appendix B in a unified notation. In the present 
calculations, we approximate C (=J r g + /*„ S-T 
1975) by including only the linear terms in F m . We 
also neglect <uj 2 >, which is relatively small over the 
time span of wave evolution we shall consider (i.e., 
5.fi days). 

Note that our mode! still contains many of the 
deficiencies that were inherent in the previous 
studies, e.g. (i) the two-level approximation elimi- 
nates the possibility for treating the true vertical 
structure of developing cyclone waves in complete 
detail, (ii) the latitudinal boundaries are “soft” in 


Tellus 37A (1985), 1 


44 


B. SALTZMAN AND CHUNG MUH TANG 


the sense that although they are closed for 
geostrophic (non-divergent) flow, they are open for 
non-geostrophic (divergent) flow (cf.. Moen, 1974; 
Paegle and MacDonald. 1974), (iii) in treating a 
single primary wave we neglect barotropic wave- 
wave. as well as wave-mean flow, interactions and, 
(iv) no feedbacks to the basic zonal thermal state 
S m ) are included as they were in S-T (1982). 
so that the primary wave grows exponentially 


without saturating (i.e., without reaching the 
occlusion stage) and we can only consider the early 
stages of wave growth during which the finite- 
amplitude wave attains its fundamental structure. 
One consequence of this latter deficiency is that, 
along with the unabated growth of the amplitude of 
the motion field, the amplitude of the mixing ratio 
field e will also tend to grow beyond realistic 
bounds. Potentially, this could even lead to nega- 





Fig. 2. Solution for fields of \y 2 and tpj in units of 10 7 m 2 s'* (left) and oj z in units of 10*° mb s'* (right), at t « 4.8, 
5.2. and 5.6 days. For t = 5.6 days, full and open dots denote closed centers of 750 mb lows and highs, respectively 
(also shown in Figs. 3, 5 and 6). 

Tellus 37A (1985), 1 


ORIGINAL PAGE IS 
OF POOR QUALITY 









FINITE-AMPLITUDE BAROCLINIC WAVES 


45 


tive values of c y For all of these reasons, our 
results can only be expected to give a qualitative 
understanding of the role of baroclinic. non -geo- 
strophic processes in the early evolution of the 
wave and its associated constituent concentrations. 


3. Solutions for the fields of passive 
constituents 

The analytical solution of (6a, b, c) (7a, b, c) and 
(8a, b) for the evolution of a baroclinic wave are 
given in S-T (1972, 1975) and will not be repeated 
here. In Fig. 2 we plot the solution, corresponding 
to the same basic state and initial conditions given 
in S-T (1972), at roughly 10 h intervals for three 
points in time after the introduction of the initial 
perturbation: t = 4.8, 5.2, and 5.6 days, the last of 
which (5.6 days) was discussed previously in S-T 
(1972. Fig. 1 la, b; 1975, Fig. 4). According to the 
solution, the wave progresses eastward at the phase 
speed c r = 16 m The analytical solutions of 
(6d). (7d), and (8c), governing the separate 
components of e y are given in Appendix C. 

Suppose we consider, as Case 1, that a non- 
interacting (passive) constituent, confined to the 
lower layer, is uniformly distributed horizontally at 
the time of inception of the wave at t = 0 and has 
an arbitrary mixing ratio, = & 3 0) = 2 x 10“ 3 . 
In Fig. 3 we show the composite distribution of c 3 
(jc,y, /) at r = 4.8, 5.2, and 5.6 days obtained by 
summing the solutions of (6d, 7d, and 8c) given in 
Appendix C, in accordance with (5). 

In Fig. 4 we show the separate components e ( 3 n , 
ey», and 4 2t that comprise the complete held shown 
in Fig. 3 for / = 5.6 days. We can see that the 
first-order eddy solution e\ l) alone would give a 
symmetrical field related to the field of divergence 
(or, equivalently, of S-T 1972) generated by 
the first order baroclinic wave. However, the 
second-order, non-geostrophic processes that be- 
come large as the wave evolves lead to a more 
asymmetric pattern in which the maximum and 
minimum concentrations are displaced to coincide 
more closely to the low and high pressure regions, 
respectively, with a relatively strong gradient 
between them roughly in the region of the thermal 
front. 

If we consider, as Case II, an initial state in 
which the constituent is not uniformly distributed, 
but is characterized by a constant equatorward 



Fig. 3. Solution for field of c 3 in units of 10~ J , for Case 1 
l* 3 °> = 2 x 10~ 3 and = 01 at / = 4.8, 5.2, and 5.6 

days. 


gradient, + 3x 10" 10 nr 1 , the first- 

order solution will contain a contribution from 
meridional advection as well as from divergence. In 
this case, the sequential solution for f = 4.8, 5.2, 
and 5.6 days is of the form shown in Fig. 5. As 
noted in Section 2, one obvious deficiency of the 
solution shown in Fig. 5 is that under the influence 
of the unabated exponential growth of the baro- 
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Fig. 4. Component solutions, e\ l \ c\ l \ and at 5.6 
days, for the complete Case I solution shown in Fig. 3. 



Fig. 5. Solution for field of e 3 in units of 10" 3 , for Case I! 
It?' - 2 X 10" 3 and dZ*'/dy *-3x 10“'® nr 1 ) at t = 
4.8, 5.2, and 5.6 days. 


clinic wave, the amplitude of e also grows at an 
unrealistic rate so that, depending on the initial 
mean value of e, negative values may even be 
generated. This occurred at f * 5.6 days. Although 
this is physically unacceptable, it is likely that if 
proper account were taken of wave occlusion 
processes (e.g~, S-T, 1982) this unrealistic feature 
would be removed without altering the general 
nature of this solution. Aside from this increase in 
amplitude, this solution for d& 0) /dy * 0 is similar in 


form and phase to the solution shown in Fig. 3 for 
d& 0) /dy « 0. 


4. Inferences concerning water vapor and 
cloud distribution in barodinic waves 

As is well-known, water vapor is not a “pas- 
sive" constituent of the atmosphere; condensation 
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and evaporation are continually taking place, and 
the condensation process entails release of latent 
heat that can modify the temperature and motion 
fields of the waves. Nonetheless, we shall assume 
that these processes are small enough relative to the 
advective processes that we can neglect them as a 
first approximation and explore the implications of 
our solutions for the distributions of water vapor 
and cloud in baroclinic waves. Thus, if we let e = c v 
(water vapor mixing ratio) we see from Fig. 5 that 
under the influence of a growing baroclinic wave, 
an initial “climatological” basic state characterized 
by QD) = 2 x 10“ 5 and dl^/By = -3 x 10’ 10 
m _> would become eddied into a high concen- 
tration pool in the low pressure region and low con- 
centration pool in the high pressure region separated 
by a strong east-west gradient in the frontal region. 
Although, for the reasons given previously, we find 
unrealistic values of the maxima and minima at 
t - 5.6 days, we assume the distributions are 
qualitatively correct, and can even imagine that the 
neglected evaporation, condensation, and preci- 
pitation processes would at least partially offset 
these unrealistic maxima and minima. That is, to 
some degree we can identify supersaturation with 
the transformation of vapor to liquid water. 

With the above limitation in mind, we next 
obtain from the mixing ratio distribution and the 
thermal field of the wave (extrapolated to level 3) a 
qualitative estimate of the held of relative humidity 
r, defined by 






(9) 


where T z = + 1 %K. For this purpose, we simply 

set e v = 0 where c* < 0 and allow r to achieve 
highly supersaturated values (> 100%) where e v is 
large. The patterns of r obtained for / = 4.8, 5.2, 
and 5.6 days are shown in Fig. 6. 

If as in many dynamical models, we assume that 
cloud amount is proportional to r (e.g., 
Smagorinsky, 1960; Vallis, 1982); we can see that 
the solution qualitatively portrays the evolution of 
the “comma-cloud” structure commonly observed 
in satellite images of developing mid-latitude 
cyclones. Further, if as in some other models 
(Lewis, 1957) we relate cloud coverage to some 
arbitrary critical values of r and also vertical 
motion (-— w 2 \ we can obtain cloud distributions 
that may exhibit even more realistic evolutions to 





Fig. 6. Field of relative humidity r, in units of %, 
corresponding to e 3 fields for Case II shown in Fig. 5, 
for/* 4.8, 5.2, and 5.6 days. 


the comma-shape. For example, in Fig. 7 we show 
the cloud coverage obtained by applying the 
conditions r > 0.60 and w < 0 to the results shown 
in Figs. 2 and 6. This sequence can be seen to 
compare favorably with the schematic cloud 
evolution sequence in a cyclone wave given by 
Anderson and Veltishchev (1973, Fig. II1.5.5), and 
with Fig. 8 showing a sequence of GOES infrared 
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condition (r > 0.60, u> <0) applied to Fig. 6, for t = 4.8, 
5.2, and 5.6 days. Dotted lines labeled A, B, and C, for 
i = 5.6 days, are those along which the vertical cross 
sections shown in Figs. 10, 1 1, and 12 are drawn. 

images spaced roughly 10 h apart (as in our 
solutions) for a cyclone-wave development over the 
eastern United States. In the last figure of this 
GOES sequence (Fig. 8d), the storm is beyond the 
occlusion stage and is strongly under the influence 
of non-adiabatic processes over the ocean. The first 
three images can be roughly compared to the 


results shown in Fig. 7 for f = 4.8, 5.2, and 5.6 
days. We note the important roles played in the 
comma development at t - 5.6 days by the 
emergence of a strong subsidence region and 
associated low-level divergence in the region just 
equatorward of the low-pressure center which 
generates the dry clear zone (see a> 2 . Fig. 2, 5.6 
days) and by the horizontal advective processes 
that sweep water vapor poleward around the 
low-pressure center at low levels. 

An excellent discussion of the synoptic features 
associated with comma-cloud development has 
been given by Carlson (1980). His analysis 
suggests that two main cloud-producing streams of 
air relative to the moving storm system are 
involved in the generation of the comma cloud. One 
is characterized as a potentially-warm “conveyor 
belt” bringing moist air from the low latitudes, in 
which middle- and high-level cloud form in the 
“tail” and eastern portion of the comma; the other 
is a potentially-coolcr low-level “conveyor belt” 
from the east in which the stratus shield con- 
stituting the western portion of the comma forms. 
In spite of the limitations of our two-level model, 
these flow features, summarized in Carlson’s Figs. 
9 and 10 (1980), seem to be well replicated by our 
solution. The upper and lower level air streams, 
relative to our eastward moving baroclinic wave 
system containing the qualitative comma cloud, are 
shown in Fig. 9 (left) along with the satellite 
infrared cloud image shown in Fig. 8c (right). Also 
shown in Fig. 9 (left) is the branching of the flow at 
level 3 (dashed streamline) near the tip of the 
comma tail, as described by Carlson (1980). In 
Figs. 10, 11, and 12 we show the vertical cross 
sections of the solution flow along the three slices 
labeled A, B, and C in Figs. 7 and 9, including a 
representation of the variations of temperature and 
moisture along section C. The low level “cold 
conveyor” is depicted in section A (Fig. 10); the 
cold front and adjacent dry-descending and warm- 
moist-ascendmg poleward air streams are depicted 
in section B (Fig. 11); the upper level “warm 
conveyor” is depicted in the north-south section C 
(Fig. 12). Although these features emerge naturally 
in detailed numerical forecasting solutions, we 
believe the present, two-layer, highly-truncated 
mode) is the “minimum” low-order system capable 
of demonstrating that these complex features 
follow, qualitatively at least, from the basic laws 
governing atmospheric behavior. 
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Fig. 8. Sequence of GOES infrared images, spaced roughly 10 h apart, of cloud coverage associated with an evolving 
cyclone-wave over the eastern United States, (a) 2105 GMT, 24 OCT, (b) 0600 GMT, 25 OCT, (c) 1500 GMT, 25 
OCT. (d) 0530 GMT, 26 OCT. 


5. Concluding remarks 

As we have noted, the present study is meant to 
complement the more detailed numerical modeling 
efforts aimed at accounting for and predicting 
large-scale weather phenomena. It is at least of 
pedagogical value to be able to show that major 
observational features can be explained in terms of 
highly idealized, simple, models that contain a 
minimum of essentia) physical ingredients. 

In the present case, we have (1) deduced the 


evolving distribution of an arbitrary low-level 
passive tracer or pollutant under the influence of a 
growing barodinic wave, and (2) applied this 
solution to show that a two-level, adiabatic, 
barodinic model including non-geostrophic, finite- 
amplitude, passive-advective effects is apparently 
suffident to yield a qualitative solution for the 
mean features and underlying dynamics of comma 
cloud formation. It is dearly possible and poten- 
tially instructive to improve upon our simple model 
by many further additions short of a full numerical 
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Fig. SO. West-east vertical cross section, (A) of Figs. 7 and 9, across the comma cloud region of the wave poleward of 
the surface low pressure center, showing the solution vector air motions relative to the eastward moving wave. Note 
the low-level “cold conveyor" upper-level “jet stream”, and “warm conveyor” streamflow (labeled “WC”) into the 
plane along cross section C of Fig. 7. The locations of the surface low and high are indicated by L and H. The comma 
cloud region is labeled on the abscissa, and a schematic conceptualization of the probable cloud cross section is shown 
by the scalloped region. 



component of the wind v (in units of 10 m sec* 1 ), the maxima of temperature T and relative humidity r. Note the 
high-level northward-moving “dry tongue”, DT, and adjacent low-level northward-moving moist warm “conveyor”, 
WC, streamflow paralleling the frontal region. The distribution of vertical motion u> at level 2 is shown in the lower 
panel. 
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Fig. 12. South-north vertical cross section (C) of Fig. 7, showing the ’‘warm conveyor’* flow overriding the east-west 
“cold conveyor” flow into the plane of the cross section marked by CC. The distributions of Tand r along this section 
are shown in the lower panel. The comma cloud extent is indicated on the abscissa, and a hypothetical cloud cross 
section is shown by the scalloped region. 


treatment. These include additional levels in the 
vertical, wave-wave and wave-meanflow inter- 
actions, wave saturating feedbacks with the basic 
state, constituent fluxes to and from the upper half 
of the atmosphere, and the replacement of our 
semi-empirical model for cloud by simple rep- 
resentations of cloud and precipitation physics 
including the non-adiabatic effects of latent heat 
release. 


6. Acknowledgments 

This study was supported by the National 
Aeronautics and Space Administration under 
contract NAS 8-34903 at Yale University and 
NAS 8-34010 at the Universities Space Research 
Association. 


7. Appendix A 

List of symbols 
x distance eastward 
y distance northward 
2 height above sea level 
p pressure 
t time 

u eastward wind speed 
v northward wind speed 
a) dp/d/ 

g acceleration of gravity 
c p specific heat at constant pressure 
R gas constant for air 
/ Coriolis parameter, 20 sin + 
fl angular speed of earth's rotation 
ft latitude 
P df/dy 
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d> geopotential, gz 
W geostrophic streamfunction,/* 1 C> 

X velocity potential 
T temperature, — pfR “ 1 dy/dp 
a specific volume 

6 potential temperature, R~ l pj p l " * a 
S static stability, — ad In 6/ dp 
k R/c p 

i unit vector eastward 
j unit vector northward 
k unit vector upward 

V i dldx +jd/dy 

V vector horizontal wind, ui + vj = V m + V x 

V 0 non-divergent (geostrophic) wind component, 
k x Vy = u 0 i + v m j (u 0 = -dy/d\\ v 0 = By// 
dx) 

V x irrotational wind component, V* = u x i + v x j 
(u x = dy/dx, V x = £*/£>’) 

f vertical component of vorticity, A*V x F = 
V 2 y 

S horizontal divergence, V - V= V 2 x 
and, with reference to Fig. 1, we also define 
Vi = + v'*). 

v't * - V',). 


V(V 2 ^ n) ) + 0) iU — V 2 
dp 

_ d<o"' ldw<" du<" 

- V 2 ur' u — 

dp \ dy dp 

dw'" dv<"\y r d d 

£x dp Jj [ £y 3y 

d 

I — —I + w (l] — V 2 Uf il) 

\ dy 2 / ty 

d (d 2 y/ a) \ d 2 y' u du> ili 

•f cu m — I ) 

dp\ df / dy 2 dp 


, to*” 

- v 2 W il) — 

dp dy dp dy dp 


r 

d0)<" dv ll) 




Pi \ By Bx 1 cx By By 1 ) 

Ifi B V ™ 8 \ 


T = 





dv T 
dy 




dx * 


8. Appendix B 

Forcing terms in eqs. ( la\ (8a) and (8b) 


9. Appendix C 

Solution for the mixing ratio e 3 
Using the same notation employed in S-T (1972, 
1975, 1982), the solutions of (6d), (7d), and (8c) 
for the components of e j(x, y, t ) represented by the 
expansion (1) can be expressed in the following 
forms, respectively: 

c5, M (x, y\ t) — (E r cos kX — E ( sin kX) sin ly • e* 

(Cl) 

t^(y\ t) m (C r cos 2 ly - G, sin 2 ly) e 2 * (C2) 


4 

A « 

Pi 

dVl 

dy 


r>v ( vV0-*-(^xi^ 

r X 3 >’ /J’ 


^ = - | F*"- v(v J ^‘")r + v 1 

£x 


0 
be 

$>■ l «>•’ ;J v * e, Sp ] 


c*, 5 ‘ = Z Z Z llY/m, n,/>) cos mly + n r (m, n,p) 

Rip 

x sin m(y] cos nit* + lY/m, n, p) cos mly 
+ n/m, n,p) sin m/y] sin nAJT } e*" (C3) 

To save space, the coefficients in (Cl), (C2) and 
(C3) can be written in terms of a generalized 
notation, whereby paired relations of the form 

Z r = —A +£,+ C, + D-fl, 

Zf = —8/4* C r + D • fc, 
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are condensed to the single relation 


where 


— +4 i ^ ^/*F ^)* 

and where the following definitions are introduced. 


£';i= 


it y o> 

k. PiM 2 V * 


(c,fi r . ( ±i/n, r ), 


1/ 

V i mV 1 * ef, 

W 1 m U 7 + 9 cl 

ij% (8f>/8y) } , 

// 1 I4' 2 I j 

A 2u'*\2SP p\+P)‘ 


As noted, all other quantities are as defined in S-T 
(1972, 1975, 1982). Thus, the coefficients in (Cl), 
(C2) and (C3) can be written in the following 
forms. 



G,= Ae<?\ 


Gj = — 2lAi">'--J-(E,'¥ il -E l 'V }r ) . 

4c ' 

The only coefficients appearing in the expansion 
(C3) for are of dimensionality (1, 1, 1), (0, 2, 2), 

(2, 2, 2), (1, 1, 3), and (3, 1, 2), that is. 


y m o, i, d=±£;* ( '± 

n,,(i, 1, 1)«±£«,», 

Y m (0,2, 2) = ±£*/±£<”». 

n r . ( (0,2,2) = 0, 

Y M (2,2,2) = ±£^' + ^£« b /, 


Il r . l (2,2.2)=±£«?', 

Y,.,(l, 1, 3) - +£«/» ± E'*?> ± £<?,«, 
Hr./d, l,3)-+£«?»±£<^"±£«/', 
Y f ( (3, 1, 3) - ±£«," ± £«,» ± £«,”, 
IKJ3, 1, 3) = ±£“?' ± £<?/> ± £«,», 


* . W,.. 

4p 2 ^ K l 

£«f/»= ±A/i«» lf'- l (3c j 'J' J( , r + IW JM ). 

£«'?» = 2Alf'- l I-3c ( (4- , y£ i>r + /* 4'j/., £?>) 

+ l/(4-/£,. ( + /> ¥*., «P)1 + M F ' *»i 

4 c t W~ 

-E.'VJQC'V^+UVi,',)], 



£«,” = _K- - V,) ej°» - (*/>,)*' c«, # > 

x (s*d, i, i), era, 1 , i»j + w(v„-j»,) e< # > 

-(W'ei # '(Cr(M. 1X5;(1. 1, 1))]), 

£««»= -K-Mc./Otpj^ 1 )"'^ 0 ' 

x (5;a, i, i),-era, 1. i»- w(*p,^)-' ?;> 
x(era,i,ixs:a, 1. i»i* 


£««> - -K-M+C^H-^j,^,) £« # ' “ (2 *?,)-' £?> 

x(Sr(2,2,2Xe:(2,2,2))l 

+ Glit l (^„ ~^ 3 ) £y W - (2kPi)~ l £<3°> 
x(e?(2, 2,2),5:(2, 2,2)))), 

£<m = -K- l I±c ( /(4Jtp J p*)-* £«•' - 

x (5J(2, 2, 2X CJ(2, 2, 2)) - W(4*p : * V e«» 
x (Or(2,2,2XS:(2,2,2))l, 
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£ ,n, m -rj £<;> + a k Pj )-' *,•> 

x (S*(Q, 2, 2), CJ(0, 2, 2))) + UluK-fr^ rj 
x *„«' + (2 kp,)-' £««' (C:(0, 2, 2), s;(0, 2, 2))J |. 

£(..»> = _H'-M+3c,|o,. ( (l, n)c«" + ( k Pl )-' 
x (<5(w) a r ( (2, n) - £«•’ a r> ,(3, «))] 

+ a,. f (4, it) + ikp } )-' (i(n) e»> a„,(5, n) 

— «j #) a,.^6, n)))},(n * 1, 2, 3, 4). 

In the above expression for £}*} n> , the quantities 
S(n) and a r J (1-6, n) are defined as follows, with 
the subscripts (r, i) deleted: 

f 3//(* 1 + 9/ J ), n=l,3 
*" )= [//(*’ + /’) n = 2, 4, 

a(l,l) = (-^3,.r,), 

a(2, 1) = l(*J(3, 1, 3) + R 9 ( 3, 1, 3)), (*T(3, 1, 3) 

+ *,(3,1,3))], 

a(3, 1) = [5*(3, 1,3), GJ(3, 1,3)], 
a(4, ])=(>T 3 ,- ^ 3 ), 

a(5, 1) = l(PT(3, 1, 3) + *,( 3, 1, 3), (*:(3, 1, 3) 

+ *,(3,1,3))], 

a(6, 1) = l£?J(3, 1,3), S;( 3, 1,3)], 
a(l, 2) = (— . *~ 3 , 3 ), 

a(2, 2) = |(*;(1, 1, 3) + *,(1, 1, 3)), (*J(1, 1, 3) 

+ *,(1,1,3))]. 


a(3, 2) = lS*(l, 1,3), £(1,1,3)]. 
a(4, 2) =(^,,-.*' 3 ), 

a(5, 2) = [(*•(!, 1, 3) + *,(1, 1, 3)), (*J(1, 1, 3) 

+ *,(1, 1,3))], 

a(6,2)=lGJ(l, 1,3),5;(1, 1,3)], 

a(l,3)=(-r„ sT,), 

a(2,3)=-a(3,l), 

a(3,3)»a(2, 1), 

a(4,3) = <r„-.r,), 

o(5, 3) = -a(6,l), 

a(6, 3) = o(5,l), 

a(l,4) = (- >r J , ^ 3 ), 
a(2, 4) = -o(3,2), 
a(3, 4) = o(2, 2), 
a(4, 4) = ( ^ 3 , - r 3 ), 
o(5, 4) = -a(6,2), 
a( 6 , 4) =a(5, 2). 

The quantities *J, Q*, *;, and SJ are equal to 
the quantities * a , g a , and S B , respectively, 
given in Table 2 of S-T (1975), minus the 
non-linear (quadratic) terms). 
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ABSTRACT 

As an extension of out previous study (Saltzman and Tang, 1985) in which 
water vapor and cloud are confined to the lower half of the atmosphere in a 
growing baroclinic wave, we now allow for the vertical flux of moisture across 
the middle (e.g., 500 mb) level to and from the upper level. The results show 
that there is a tendency for the upper-level maximum mixing ratio to shift to 
the east relative to the lower-level maximum mixing ratio. This is consistent 
with the eastward tilt with height of the common form cloud in the “cold- 
conveyor” cross-section and the poleward tilt of the cloud in the “warm- 
conveyor” cross-section as observed by Carlson (1980). 


1. Introduction 

Saltzman and Tang (1985) (hereafter 
referred to as ST— 85) considered the 
effect of finite-amplitude baroclinic 
waves on passive, low-level, atmospheric 
constituents. In particular, if water 
vapor is considered to be such a passive 
constituent it was shown that the 
implied relative humidity and cloud 
distribution would tend to evolve into 


the comma form commonly observed 
on satellite images of mid-latitude 
cyclone wave. The wind, temperature, 
relative humidity, vertical motion, and 
comma cloud fields exhibited in the 
cross-sections seem to agree with the 
salient features (e.g., ‘‘cold conveyor”, 
“dry tongue”, and “warm conveyor”) 
deduced from observations by Carlson 
(1980). Since the moisture field in the 
model was confined to the lower layer 
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only, the schematic cloud boundaries 
above the middle level was speculative 
(cf. ST-85, Figs. 10 and 12). The 
purpose of this note is to show that tliis 
feature is reflected in the upper-level 
mixing ratio field in the model when the 
upper-level moisture field is taken into 
account. 

2. The Model 

We consider a two-level model as in 
ST— 85, but now in addition to the 
lower-level mixing ratio field, we 
consider variations in the upper-level 
mixing ratio field. The definitions of 
the symbols are listed in Appendix A, 
ST— 85. The continuity equation for an 
atmospheric constituent is given by 

de d 9 co 

+ V • Ve + — (ecu) — e — =S (1) 

at ap dp 

Applying this equation at levels 1 
and 3 , we obtain 


+ ^(e a -e,)=S, (2a) 

dt pj 


?li+V, - Ve, + — (e,-e a )-S, (2b) 

at Pi 

Next, we assume that the mixing ratio 
at the middle level is given by 


e a « (e, — e,)/2 (3) 

which .would be the case, for example, if 
• has a parabolic profile given by 


K 

€ = — P J (4) 

Pj 

where K is a function of x, y, t, only. 
Substitution of (3) in (2a, b) leads to 

— + V,.Ve, +^-(€ 3 - 36 , ) — =S,(5a) 
dt 2 Pj 

— — + Vj • V + ^(e 3 +c.)-^ = S 3 (5b) 

at 2 p a 

Next, as in ST— 85, we (i) assume 
Si ,3 to be zero (or, alternately, that e 
represents total water content, vapor 
plus cloud), (ii) expand each of the 
dependent variables £ = (ip, x, co, e) 
into a zonal mean part £(y,t) and an 
eddy part £ ', i.e., 

£= £ + £', 

and (iii) resolve each of these parts into 
an initial, basic zonal field denoted by 
the superscript ( 0 ), a primary, quasi- 
geostrophic field denoted by the super- 
script (1) (Phillips, 1954), and a second- 
order, non-quasigeo st rophic field 

denoted by the superscript ( 2 ); i.e., 

I, l <0, + {(!)+ |(2) 

r. *«>+{«*. 


The systems of equations are as 
follows: 




(i) Equations for the primary wave 
a) ,cu (, >, x (1) ,e (I) J: 
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S W 2 fli 2f «>„. W 
- V cj' 1 ^ w, + X- 

f 2 2 2 P> 9x 

P 2 


( 1 ) _<i) 

a a Xi, 3 _ - 5, 

3y J p 2 


, 4 _(0) 

3 


+ — u T 

P2 

dx 

v 2 ^ 



a) 

„ 2 (1) 


« 2 

V *1.3 - 

+ - 

* 

P 2 


(6b) Be/ 1 * 3 (i) (i) (i)3e (0) 

= (v. e ), - (u v X 

3t 3y $ x 3y * 

(6c) - - — (e 3 - 3e, ) Wj , 

2Pi 


* (1) * «> 
3e, _©) 3e, 

— +u. 

3t 3x 


3\f'{ 1> dc{ 0) 1 ..(0) ,_(0> 

— — (e 3 — 3e, ) 

3x dy 2p 2 


a =0) a a^O) 

1 «®>4. =©) v (1) 

(e 3 +e, )cjj 

2p 2 


(III) Equations for the secondary 
wave fields [ ^ (2) , w (2) , e (2) J : 


(l) CD 

de 3 (0) 3e 3 

3t + Ua d 

._a^ 

3x 3y 


1 ,_©> _©K (1) 

t — ( e 3 +c, )w 2 

2pj 




(II) Equations for the forced zonal 
state ($ u \ c«j w \ x (1 >, e* 1 *]: 

S jtoP> 2f _ m m 

r^^" fA( * k < 7 ’> 


3 3 J Ml) 

— ( — — ) 

3t dy* U 

•-ivj*. V(V** n )] 1> j±^C I “ > (7b) 


Pj 

(8a) 


Pa 

(8b) 


A £ p» +( rfo.^, 

at ax 


a2 ©) a,.UJ a_UJ 

+( u;>!L ) +( ,j >, ‘ )■ 

* by A 3x * ay 


a*«> a««> *=«) 

+ u({)il ) + v a)»f + «»»£_ 
dx v dy x 3y 
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+ v 


0)8f®+ 
^ dy 


aJ (0) 

4 2) j 

X ay 


2Pi 

+ (80 


dt 


£, +(u — >3 

dx 

— t V ‘ 1 ’ — ♦ <u , '> £ ) •+ CV 

x ay x ax x a x 

.,„3e<>' .Dii" 

*“*Tx v * 57 x ay 


A <2) 

+ v ^ 


3e (0 > 


( 2 ) ( 2 ) 

* ay 


} 3 


2p 2 

(8d) 


As in ST— 85, where the double 
signs appear the upper sign refers to level 
1 (250 mb) and the lower sign refers to 
level 3 (750 mb); $ is the horizontally 
uniform value of static stability, and A, 
B, and C are the non-homogeneous 
forcing functions due to lower order 
fields given in Appendix B in ST-85. 
Again, we approximate C (=Fg + F a , 
Saltzman and Tang 1975) by including 
only the linear terms in Fg. We also 
neglect u> 2 (2). The solutions forei and 
e 3 can be obtained by using the formulas 
for the streamfunction, temperature, 
vertical motion, and velocity potential 


in Saltzman and Tang (1972, 1975), 
plus the basic mixing ratio and its 
meridional gradient. The derivations 
are straightforward, but rather lengthy, 
and we shall not present them here. 

3. Results 

Before we present the new solution, 
we show in Fig. 1 a corrected version of 
Fig. 5 of ST-85 portraying the 
evolution of the low-level water vapor 
mixing ratio field ej . Notice that the 
main error is in the region of sinking 
motion where a weaker minimum 
appears, corresponding to only a small 
area of negative values. Since the 
western cloud boundary is mainly 
determined by the eastern boundary of 
the sinking region, the comma-type 
cloud distribution remains virtually the 
same. 

The parameters in the present model 
are the same as those in ST— 85, except 
that we consider only the case <F, (0) 
(D/2)=(2/9)xl0-*, e, (0) (D/2)=2xl0-» 
and de, (0) /3y ® — (1 /3)x 10* 10 m-» , 
ae 3 (0, /dy®— JxlO'^m* 1 . Note that 
these prescribed values satisfy (4). 

In Fig. 2 we show the mixing ratio 
at level 1 (250 mb) and at level 3 (750 
mb) at t*4.8, 5.2 and 5.6 days. The 
upper-level maximum tends to be shifted 
eastward relative to the lower-level 
maximum as would be expected given 
the increasing westerly wind speed with 
height. The lower-level field is similar 
to Fig. 1 but now e3 is everywhere 
positive and its maximum is smaller than 
that in Fig. 1 . However, negative values 
appear at the upper level. As noted by 
Lorenz (1984) this feature arises from 
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here we simply accept only the qualita- 
tive result that these negative values 
indicate dry regions of vanishingly small 
water vapor content. 

In Fig. 3 we show the relative 
humidity field at the lower level 
estimated qualitatively by 




Fig. 1 Solution for field of e 3 in units of 
10-3, f OT e ,=0. e,(°)=2xl0-3, and 
Sej(°)/dy* —3xl0-10,n-l f at t*4.8, 
5,2, and 5.6 days. 


the fact that the horizontal variations of 
€ are comparable to its mean value and 
cannot be represented realistically by 
only a few spatial components. In lieu 
of adopting Lorenz’s more rigorous 
suggestions for avoiding this difficulty, 
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Fig. 2 Solution for the fields of £j (left) and 
? 3 (right) in units of 10*3, for £, (0)= 
(2/91x10-3, e-(0) . 2x10*3, and 
di t (y)/dy * -<l/3)x lO^m-l, and 
d£j(0)/3y=— SxlO' 10 ™*!, at t»4.8, 
5.2, and 5.6 days. 


We note that unrealistic supersaturation 
still occurs. In Fig. 4 we show the cloud 
coverage obtained by applying the 
condition used previously, i.e., r> 0.50 
and 03 < 0. As before the comma or 
scallop shape cloud distribution is seen 


to develop as time progresses. 

4. Conclusion 

It appears that the qualitative cloud 
sketch in the cross-sections shown in 
Figs. 10 and 12 of ST-85 is supported 


r 

4M 




Fig. 3 Field of relative humidity at the lower 
level, r, in units of %, corresponding 
to £ 3 fields shown m Fig 2, for 
t*4.8, 5.2 and 5.6 days. 
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Fig. 4 Cloud areas deduced from the condi- 
tion (r ) 0.50, U) < 0) applied to Fig. 3, 
for t»4.8, 5.2 and 5.6 days. 


remain however. 
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by the present calculation of the 
moisture fields in the upper and lower 
levels. Note,, in particular, that the 
eastward tilt of the cloud in Fig. 10 of 
ST— 85 is consistent with the eastward 
tilt of the maximum mixing ratio in 
Fig. 2 of the present paper. The other 
deficiencies mentioned in ST-8S still 
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ABSTRACT 

Data from the National Meteorological Center and net outgoing infrared radiation (IR) data measured by 
NOAA satellites for January 1977 are used to compute estimates of the spectral and spatial contributions to 
the net generation of available potential energy in the Northern Hemisphere due to infrared radiation. Although 
these estimates are necessarily crude, the results obtained indicate that IR causes destruction of both zonal 
and eddy available potential energy. The contributions from midlatitudes to the zonal and eddy generation 
are Gz m —5.0 W m' J and G c •* —0.6 W m* J , respectively. The eddy generation is due almost entirely to 
stationary wavenumbers 1 and 2. Comparison with earlier studies and computation of Newtonian coaling 
coefficients are discussed. 


1. Introduction and background 

Diabatic heating in the atmosphere is composed of 
four principal contributions. These are: 1) absorption 
of solar radiation; 2) emission and absorption of long- 
wave (infrared) radiation; 3) heating due to the release 
of latent heat of condensation; and 4) transfer of sen- 
sible heat from the surface of the earth to the atmo- 
sphere. The characteristics of atmospheric diabatic 
heating are of crucial importance because it is the non- 
uniform distribution of diabatic heating that is the 
driving force behind the general circulation. A measure 
of this effect in the atmospheric energy cycle is given 
by the generation of available potential energy (APE), 
which is determined by the spatial correlation of the 
diabatic heating with the atmospheric temperature field 
(Lorenz, 1955). The rate of APE generation is one 
measure of the intensity of the general circulation 
(Wiin-Nielsen, 1968). 

Diabatic heating and APE generation are difficult 
quantities to measure and are conventionally estimated 
as budget residuals. Two approaches have been used 
traditionally to compute the diabatic heating. First, 
the heating is computed as a residual in a heat budget 
(Wiin-Nielsen and Brown, 1962; Brown, 1964; Geller 
and Avery, 1978) or potential vorticity budget (Lawn- 
iczak, 1969). In principle, the residual methods can 
give the spatial and temporal distributions of the beat- 
ing field. However, they are limited by data quality 
and the need to compute vertical velocities in the heat 
budget calculation. The second approach is to use de- 
tailed radiation models to compute the time mean and 


1 Present affiliation: Meteorology Research Center, Control Data 
Corporation, Minneapolis, MN 55440. 


(often) zonal mean radiational heating fields (Katay- 
ama, 1967; Newell et a!., 1974; Freeman and Liou, 
1979). Latent beat release and surface sensible heat 
flux terms must be estimated independently. Unfor- 
tunately, this approach must necessarily deal with the 
time mean heating because the radiation models re- 
quire assumptions about the mean temperature, 
cloudiness (amount and level) and atmospheric com- 
position (water vapor, aerosols, etc.). Therefore, these 
studies by their very nature cannot deal with the tran- 
sient part of the diabatic heating and APE generation. 
Newell et al. (1974) were further constrained to the 
time averaged, zonal mean APE generation, although 
Katayama (1967) made estimates of the standing eddy 
APE generation. 

In this report, we use an alternative approach to 
estimate the contribution of infrared cooling to the 
APE generation. We use NOAA satellite measurements 
of outgoing infrared radiation (Gruber and Winston, 
1978) combined with NMC height field data to de- 
termine temperature, both for January 1977, to esti- 
mate the spectral and geographic contributions to this 
component of the generation. Suomi and Shen (1963) 
were the first to use this approach to estimate IR- 
induced APE generation on selected dates and for se- 
lected areas dtfring 1959 and 1960. They found the 
generation due to IR to be comparable in magnitude 
with other estimates of the total APE generation and 
to exhibit significant temporal variability. 

2. Data and approach 

The total rate of APE generation Gcan be expressed 
as the sum of various spatial and temporal contribu- 
tions: 

G * Gz + Gg * Gzs + Gzt + Ges + Get, 
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where subscripts Z refer to zonally averaged quantities 

1 C 2w 

< < * 

and E to eddy quantities ( ) E = ( ) — ( ) z , and 
the subscripts 5 and T refer to the stationary (time 
mean) mode 

1 f* 

1 * 

and transient mode ( ) r = ( ) - ( ) s , respectively. 
Using Lorenz’ (1953) approximate form, these con- 
tributions can be written as 


Gzs “ VlQzsT'zs\, 

(1) 

Gzt m [7(9 'ztT'zi)s[, 

(2) 

Ggs “ MtoT'sslz]. 

(3) 

Get m [7 (QetTbt)zs]> 

(4) 


where the brackets denote an average over the mass 
of the atmosphere normalized per unit area, q is the 
diabatic heating rate per unit mass, T the temperature, 
7 “ tf/ifcpc), 
and 

_ _ add 

° m 6 dp 

is a measure of the static stability. A value of a for 
the 850-300 mb layer was adapted from Tomatsu 
(1979). Departures of the zonal mean heating and 
temperature from their hemispheric mean values are 
denoted by q'z and 7z, respectively. (We will ignore 
departures of the hemispheric means from the global 
means.) Dutton and Johnson (1967) and Newell el al. 
(1974) discuss the validity of these approximate forms 
of Gzs and Gzr • In addition, the eddy components of 
the generation can be broken down into zonal har- 
monic spectra 

Gg m 2 “ Z [ yiQmTm 4" 9m 

m m 

where the subscript m denotes the contribution for 
wavenumber m, and an asterisk denotes a complex 
conjugate of the coefficients q m or T m . Thus, relative 
cooling of relatively warm air or relative heating of 
relatively cold air will result in destruction of APE and 
vice versa. 

The infrared radiation data that we use are those 
measured by the scanning radiometers aboard NOAA 
operational polar-orbiting satellites. Descriptions of the 
instruments and data reduction procedures are given 
by Gruber and Winston (1978). The satellite’s equator- 
crossings occur at 0900 and 2100 LST daily. We av- 
erage the twice-daily observations to get daily average 


values at each grid point of a 2.5° latitude-longitude 
array. The NOAA satellite observes IR in the atmo- 
spheric “window region” (10.5-12.5 wavelengths). 
In the analysis procedure used by NOAA, estimates 
are made of the total outgoing IR flux from these 
window region radiance measurements by using a 
nonlinear regression model derived from radiation 
calculations for 99 different atmospheres covering a 
wide variety of temperature and moisture conditions 
as well as clear and cloudy skies (Gruber and Winston, 
1978). 

Let us briefly discuss the sources of the net outgoing 
IR. The surface of the earth or of a cloud radiates 
essentially as a blackbody. A portion of this radiation 
is transmitted by the atmosphere to space and the 
remainder is absorbed by atmospheric gases, primarily 
water vapor, carbon dioxide and ozone, or by clouds. 
These constituents then reradiate IR at an intensity 
appropriate to their temperature and emissivity. For 
a long-term global mean, most of the IR leaving the 
top of the atmosphere comes from the atmosphere 
itself. However, the distribution of the total outgoing 
IR is particularly dependent upon the distribution of 
cloudiness and water vapor in the atmosphere (Wein- 
stein and Suomi, 1961; Paltridge and Platt, 1976; Free- 
man and Liou, 1979; Ohring and Gruber, 1983). 

. Our interest is in estimating the tropospheric IR 
cooling rate from the satellite measurements. Earlier 
studies (e.g., Kuhn and Suomi, 1960; Sabatini and 
Suomi, 1962) indicate that satellite measurements of 
outgoing IR give a good estimate of the net IR flux 
divergence from the troposphere compared to direct 
radiometersonde measurements. However, some un- 
known contribution to the satellite-derived IR flux will 
be provided by emission from the earth’s surface. For 
a long-term, global mean this fraction will be roughly 
20% of the total (Paltridge and Platt, 1976). In cloud- 
free conditions it will be closer to 50% (or greater in 
regions of low water vapor content, i.e., near the pole), 
while in regions of overcast it will be near zero. In 
winter, average daytime cloudiness typically exceeds 
75% over midlatitude oceans and 50% over midlatitude 
continents (except in desert regions) (Clapp, 1964), 
suggesting a generally smaller contribution from the 
surface to the satellite IR observations. However, the 
vertical distribution of the cloudiness is crucial in de- 
termining its effect on the tropospheric cooling rate 
for a given temperature and moisture profile (Cox, 
1969). Unfortunately, the available satellite observa- 
tions cannot determine the vertical structure of the 
cooling rate. As an approximation, we shall assume 
that a linear relationship exists between measured out- 
going IR and the resultant tropospheric cooling rate 
due to this radiation. 

Following Suomi and Shen (1963), we will assume 
that the satellite-observed IR is representative of the 
vertical mean value for the troposphere ( 1 000- 1 00 mb 
layer) and will ignore correlations between q and T 
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Fic. I. The amplitude of [IflJ], for the total IR (solid line) and 
time-mean contribution ([l« TO l], dotted line) for January 1977 as the 
function of the zonal wavenumber m. In this case, q m is averaged 
over the zone from 30 to 80°N. 
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due to variations about their vertical average. This is 
a crude but acceptable first approximation. Freeman 
and Liou’s (1979) model, for example, indicates that 
variations in q u about its vertical mean are not ex- 
cessively large for January conditions, particularly in 
extratropical latitudes. By making this assumption, we 
will overestimate the tropospheric cooling rate by ig- 
noring the contribution to the satellite-observed IR 
flux given by radiation from the earth’s surface trans- 
mitted by the atmosphere and by radiation from layers 
above the troposphere. The former contribution could 
account for up to 50% of the total observed irradiance 
in certain locations and the latter effect is probably 
small. 

Tropospheric mean temperatures are computed hy- 
drostatically from the NMC operational analyses of 
the 1000 and 100 mb geopotential heights. These data 
are available twice daily at 0000 and 1200 GMT and 
were averaged to give daily mean temperatures. Certain 
inaccuracies may be introduced because the satellite 
data and temperature data are not synoptic, leading 
to misinterpretation of the transient part of G. How- 
ever, for slowly evolving, large spatial scales, this prob- 
lem should not be serious. 
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Fig. 3. The integrands of G z (a) and G £ (b) as 
functions of latitude. 


Although our admittedly crude assumptions may 
introduce certain systematic errors into our calculation, 
if at least the phase of the net cooling rate is correct, 
we can proceed with our estimate of the APE gener- 
ation. However, we must keep in mind that our es- 
timates of G will at best achieve the correct sign and 
the correct order of magnitude. 


Table I. Avenge values during January 1 977 for the contributions 
to the IR-produced generation of available potential energy (W m" J ) 
in the zones noted. For comparison, Katayama’s (1967) values from 
a radiation model calculation are included. 


Zone 

G z 

Gzs 

Gz 

Gzs 

30-80°N 

—4.98 

-4.96 

-0.61 

-0.61 

0-90’N 

-4.74 

-4.73 

-0.33 

-0.31 

JUtaymms (1967) IR 

(0-90°N) 

— 

-3.443 

— 

—0.242 

JCataymma (1967) SR* 

(0-90* N) 

— 

2.914 

— 

-0.005 


* Solar radiation. 


3. Results 

We will present results for January 1977 only. Cal- 
culations for January 1978 showed qualitatively similar 
features. First, consider the amplitude spectrum of the 
eddy infrared cooling rate [l^mlL (Fig. 1). Notice that 
maximum values are obtained at planetary-scale 
wavenumbers, suggesting the important sensitivity of 
q m to the large-scale temperature field. Examination 
of the geographic distribution of q £S (Fig. 2) shows 
that certain geographic features leave a strong signature 
in the radiation field. Features such as the Tibetan 




Fig. 4. The spectrum of G £ for (a) our results for the IR contribution 
in January 1977 where the solid line is the totaj G m and the dashed 
line is Gm,, and (b) Brown's (1964) estimates of the total G m averaged 
for three Januaries (1959, 1962, 1963.) Our values are averaged for 
30-80’N and Brown's values are for 20-87.5 °N. 
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plateau and Greenland are clearly evident. Most of 
these features are represented by wavenumbers higher 
than 3, however. At the planetary scale, the IR dis- 
tribution mirrors the temperature field, as we will show 
shortly. The monthly mean field of q ES for January 
1977 does show certain large departures from the mean 
for the eight Januaries from 1975 through 1982 (not 
shown). In particular, the strong cooling off the Cal- 
ifornia coast in January 1977 has an amplitude nearly 
twice the eight-January mean value. 

In Fig. 3 the integrands of G z and G E are presented 
as functions of latitude. The monthly mean values are 
contributed almost entirely by the stationary mode 
(See Table 1). The relative maximum in the G z in- 
tegrand (Fig. 3a) near the equator is due to persistent 
high cloudiness near the ITCZ and the resultant re- 
duction in radiation from the troposphere (during 
summer, preliminary estimates indicate that this effect 
can result in a positive contribution to G z in the tropics 
due to IR). A minimum in the integrand of G z occurs 
at ~12.5°N due to relatively cloud-free conditions 
and large IR cooling rates in this zone (Gruber and 
Winston, 1978). From a zero value near 30 e N, the G z 
integrand decreases monotonically toward polar re- 
gions. Contributions to G E (Fig. 3b) are only important 
between 35 and 75°N, reaching values of approxi- 
mately - 1 W m -2 between 45 and 60°N. Katayama 
(1967) estimated Gzs and G& for the Northern Hemi- 


sphere based upon mean heating rates from a detailed 
radiation model. His results are quite similar to ours 
(Table 1), suggesting that the approximations we have 
made are acceptable. In addition, Katayama's results 
for Gss due to absorption of solar radiation indicate 
that this contribution is extremely small (<0.01 W 
m -2 ). Thus, the infrared contribution to G ES may be 
equivalent to the total radiational contribution. 

The spectrum of G E (Fig. 4a) reveals that the total 
G e is contributed almost entirely by the standing 
component of wavenumbers 1 and 2. This contribution 
by the ultralong waves is competitive with the total 
G m at wavenumbers 1 and 2 computed by Brown 
(1964; January values reproduced in Fig. 4b). The small 
positive values of G m at intermediate wavenumbers 
may be the result of a cumulative effect due to positive 
APE generation in the developing stage of cyclones 
(Vincent et al., 1977), but it is too small to be of 
importance to the mean energetics. We should add, 
of course, that the spectrum of G E is not linearly related 
to the spectrum of the APE (not shown). 

The spectral energetics for the winter of 1976-77 
have been presented by Chen (1982). Comparison with 
his calculations of the zonal tp eddy APE conversion 
C(A Z , A e ), the eddy APE to eddy kinetic energy con- 
version C[A e , K e ), and wave-wave interactions for the 
ultralong waves shows that the IR-induced APE de- 
struction is a significant term in the APE budget of 
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the standing ultralong waves, accounting for 70-80% 
of the required balancing of the APE budget residual 
at these length scales. 

Because the standing ultralong waves make the ma- 
jor contribution to the IR-induced G £ , it is of interest 
to examine the geographic distribution of the ultralong- 
wave 1R heating and temperature fields, as well as their 
correlation (Fig. 5). Notice the tendency for the IR 
heating and temperature fields to be out of phase, lead- 
ing to large negative contributions to G £S over interior 
Asia with lesser contributions over the northeast Pacific 
and over western Europe and the eastern Atlantic 
(Fig. 5c). 

The destruction of eddy APE in midlatitudes mea- 
sured by Brown ( 1 964) may be largely due to air mass 
modification over the Gulf Stream and Kuroshio. This 
process would presumably take place at primarily in- 
termediate scales. If Brown’s and our results are rep- 
resentative and comparable, we could infer that a sig- 
nificant amount of the total destruction of eddy APE 
at wavenumbers 1 and 2 is due to radiational processes, 
while radiational processes are energetically unimpor- 
tant at intermediate scales. 

Although the transient mode makes no contribution 
to the total monthly mean G z and G £ , this does not 
mean that it is zero on any given day. The time series 
for G z and Gg averaged over 30-80°N (Figs. 6a and 
b) show significant departures from the monthly mean 
that appear systematic rather than random. The de- 
partures of Gg and Gg from their monthly mean values 
are nearly always directly out of phase with the tem- 
poral variations of A z and Ag (not shown), a not al- 
together surprising result which is consistent with the 
strong dependence of the cooling rate on temperature. 

One source of systematic error in our calculation is 
probably the lack of vertical resolution in our com- 
putations. For instance, we have compared the APE 
calculated from the 1000-100 mb thickness with that 
calculated from the NMC analyses of observed tem- 
peratures at the ten mandatory levels in the tropo- 
sphere. Although the APE calculated in the present 
study is fairly accurate for wavenumbers 1, 2 and 3, 
it is consistently underestimated by a factor of 3 or 
more for wavenumbers 4 and higher. A phase reversal 
of the intermediate-scale temperature waves between 
the lower and upper tr o posphere most likely explains 
this discrepancy. This implies that G m for these wave- 
numbers is also underestimated. Although the inter- 
mediate scale G m would then assume larger positive 
values, it would still be small compared to other APE 
budget terms. 

Finally, we can use our results to compute a New- 
tonian cooling coefficient a in a manner analogous to 
that used by Saltzman (1973). This coefficient, which 
is a time constant for the destruction of eddy tem- 
perature variance due to longwave radiation, can be 
defined as 


(QeTe) is 

“ [C/T E 2 ) U ] * 

We get a se 3 X 10 -7 s -1 (a -1 38 days). 

As already mentioned, Katayama (1967) indicates 
that the contribution of solar radiation to G £S is 
negligible and so its contribution to the total Gg is 
probably small also. Our a may therefore be repre- 
sentative of the total radiational damping coefficient 
and would apply only to the ultralong waves, since 
the radiational contribution to Gg for smaller scales 
is negligible. Wiin-Nielsen et al. (1967) obtained a 
somewhat larger value for a (4 X 10" 7 s" 1 ) based on 
Brown's (1964) total diabatic beating. Saltzman (1973) 
obtained a much larger value, « * 2 X 10 _t s"\ ap- 
propriate to the total diabatic beating at 850 rob. This 
seems plausible due to the large influence of boundary- 
layer effects at that level. 




Jonuory 1977 

Fig. 6. The time tenet of the IR-produeed G z (a) and G z (b). The 
values shown are the contributions to G z and G z for the zone 30- 
80°N. The dotted line indicates the value of G» and G& in (a) and 
(b), respectively. 
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4 . Conclusions 

The infrared cooling contribution to G z and G E in 
winter based on satellite observations of the infrared 
cooling rate is negative. Further, the satellite-based es- 
timates compare well with available radiation-model- 
based estimates despite the assumptions about the 
source and vertical uniformity of the observed cooling 
rate. The contribution of G z in winter is negative at 
all Northern Hemisphere latitudes while contributions 
to G e are only significant in the 35-70°N zone and 
are due almost entirely to standing wavenumbers 1 
and 2. A Newtonian cooling coefficient of -^3 X 10~ 7 
s~' applicable to the ultralong waves was calculated, 
corresponding to a radiative damping time scale of 
—38 days. A similar value was obtained for January 
in both 1977 and 1978. 

Based upon a comparison of our results with those 
of Wiin-Nielsen and Brown (1962) and Brown (1964), 
it appears that the contribution to the destruction of 
eddy APE in midlatitudes during winter by radiational 
processes is significant for the ultralong waves, while 
destruction of eddy APE by heat flux from the earth’s 
surface may be predominant at higher wavenumbers 
(m > 3). The role of latent heat release in midlatitude 
G e is still an open question (Hayashi and Golder, 1981). 

Acknowledgments. We would like to thank Professor 
Barry Saltzman and Professor Ronald Smith for useful 
discussions. This study was supported by the National 
Aeronautics and Space Administration under Grant 
NAS 8-34903 at Yale University. 

REFERENCES 

Brown, J. A., 1 964: A diagnostic study of tropospheric diabatic beating 
and the generation of available potential energy. Tellus, 16, 
371-388. 

Chen, T.-C, 1982: A further study of spectra) energetics in the winter 
atmosphere. Mon Wea. Rev., 110, 947-961. 

Clapp, P. F., 1964: Global cloud cover for seasons using Tiros 
nephanalysis. Mon. Wea. Rev., 92, 495-507. 

Cox, S. K., 1969: Radiation models of midlatitude synoptic features. 
Mon. Wea. Rev . 97, 637-651. 

Dutton, J. A., and D. R. Johnson, 1967: The theory of available 
potential energy and a variational approach to atmospheric 
energetics. Advances in Geophysics. Vol. 12, Academic Press, 
333-436. 

Freeman, K. P., and K.-N. Liou, 1979: Climatic effects of drrus 
clouds. Advances in Geophysics, Vol. 21, Academic Press, 23 1 - 
287. 


Geller, M. A., and S. K, Avery, 1978: Northern Hemisphere distri- 
butions of diabatic beating in the troposphere derived from 
general circulation data. Mon. Wea. Rev., 106, 629-636. 

Gruber, A., and J. S. Winston, 1978: Earth-atmosphere radiative 
beating based on NOAA scanning radiometer measurements. 
Bull. Amer. Meteor. Soc., 59, 1 570- J 573. 

Hayashi, Y., and I. Golder, 1981: The effects of condensational 
beating on mid-latitude transient waves in their mature stage: 
Control experiments with a GFDL general circulation model. 
/ Atmos. Set.. 38, 2532-2539. 

Katayama, A., 1967: On the radiation budget of the troposphere 
over the Northern Hemisphere (III) — Zonal cross-sections and 
energy consideration. J. Meteor. Soc. Japan , 45, 26-38. 

Kuhn, P. M., and V. E. Suomi, I960: Infrared radiometer soundings 
on a synoptic scale. J. Geophys. Res., 65, 3669-3677. 

Lawniczak, G. E., Jr., 1969: On a multi-layer analysis of atmospheric 
diabatic processes and the generation of available potential en- 
ergy. Pta.D. dissertation. University of Michigan, Ann Arbor, 
Ml pp. 

Lorenz, E. N., 1955: Available potential energy and the maintenance 
of the general circulation. Tellus. 7, 157-167. 

Newell, R. E., J. W. Kidson, D. G. Vincent and G. J. Boer, 1974: 
The General Circulation of the Tropical Atmosphere and In- 
teractions with Extra-Tropical Latitudes, Vol. 2. The MIT Press, 
371 pp. 

Ohring, G. t and A. Gruber, 1983: Satellite radiation observations 
and climate theory. Advances in Geophysics, Vol. 25, Academic 
Press, 237-304. 

Paltridge, G. W,, and C. M. R. Platt, 1976: Radiative Processes in 
Meteorology and Climatology. Elsevier, 318 pp. 

Sabatini, R. R., and V. E. Suomi, 1962: On the possibility of at- 
mospheric infra-red cooling estimates from satellite observations. 
J. Atmos. ScL 19, 349-350. 

Saltzman, B., 1973: Parameterization of hemispheric heating and 
temperature variance fields in the lower troposphere. Pure Appi 
Geophys., 105, 890-899. 

Suomi, V. E., and W. C, Shen, 1963: Horizontal variation of infrared 
cooling and the generation of eddy available potential energy. 
J. Atmos. Sci., 20, 62-65. 

Tomatsu, K., 1979: Spectral energetics of the troposphere and lower 
stratosphere. Advances in Geophysics, Vol. 21, Academic Press, 
289-405. 

Vincent, D. G., G. B. Pant and H. J. Edmon, Jr., 1977: Generation 
of available potential energy of an extratropical cyclone system. 
Mon. Wea. Rev., 105, 1252-1265. 

Weinstein, M., and V. E. Suomi, 1961: Analysis of satellite infrared 
radiation measurement on a synoptic scale. Mon. Wea. Rev., 
89, 419-428. 

Wiin-Nielsen, A., 1968: On the intensity of the general circulation 
of the atmosphere. Rev. Geophys., 6, 559-579. 

, and J. A. Brown, 1962: On diagnostic computations of at- 
mospheric beat sources and sinks and the generation of available 
potential energy. Proc. Ini. Symp. on Numerical Weather Pre- 
diction , Tokyo, Meteor. Soc. Japan, 593-613. 

, A. Vernekar and C. H. Yang, 1967: On the development of 

barodinic waves influenced by friction and beating. Pure Appi 
Geophys., 68, 131-161. 



Tellus ( 1 984), 36 A, 52-63 


A comparison of the spectral energy and enstrophy 
budgets of blocking versus nonblocking periods 


By ANTHONY R. HANSEN, Department of Geology and Geophysics , Yale University, P.O. Box 6666 , 
New Haven , Connecticut 06511 , U.S.A. and ALFONSO SUTERA, Center for the Environment and 
Man, Inc., Hartford, Connecticut 06230, U.S.A. 


(Manuscript received October 18, 1982; in final form March 7, 1983) 


ABSTRACT 

The time mean spectral energy and enstrophy budgets for blocking versus nonblocking periods 
from the winters of 1978-79 and 1976-77 are calculated. The major differences in the 1978-79 
cases were in the nonlinear interaction terms. A pronounced upscale cascade of kinetic energy 
and enstrophy from intermediate to planetary-scale wavenumbers during blocking was found. 
During December 1976, the energy and enstrophy budgets of a case of Atlantic blocking were 
very similar to the cases in 1 978-79, but quite different from the persistent, greatly amplified, 
planetary-wave pattern that followed in January and February 1977 (as identified by Chamey et 
al., 1981). Planetary-scale barodinic processes were greatly increased during the 1977 event 
along with reduced intermediate-scale baroclinic activity and the absence of the upscale 
enstrophy cascade noted in the other blocking cases. 

The predictability time (Lilly, 1970) based on the enstrophy flux function showed an increased 
predictability for the January-February 1977 event but no significant differences between the 
other blocking cases compared to the nonblocking sample. However, the reversal of the low 
wavenumber enstrophy cascade during blocking does suggest that blocking may be more 
persistent (due to reduced dissipation of the large-scale circulation) and therefore more 
predictable. 


1. Introduction 

The term “blocking” has been used to describe a 
variety of circulation features with characteristic 
time scales larger than the synoptic time scale. The 
importance of blocking events in causing extremes 
in local surface weather and the expectation that 
the existence of these persistent features would 
allow greater forecasting skill in numerical weather 
prediction models has stimulated a great deal of 
interest in them in recent years. 

Diagnostic calculations have usually identified 
blocking with large-scale baroclinic processes (e.g., 
Paulin, 1970; Murakami and Tomatsu, 1965). 
Recently, Hansen and Chen (1982) were able to 
identify the initiation of two cases of blocking with 
enhanced baroclinic energy conversions. In one 


case, intense baroclinic cyclone waves forced the 
development of a blocking high through nonlinear 
interactions. In the other case, large-scale baro- 
dinic processes accompanied the block's develop- 
ment. 

In the present study, we will attempt to identify 
any systematic differences in the time mean 
spectra] energy and enstrophy budgets of blocking 
events compared to nonblocking periods using data 
from two recent winters (1976-77 and 1978-79). 
Recently, Bengtsson (1981) has shown from 
numerical simulations that during a blocking event 
the atmosphere is in general more predictable than 
during nonblocking periods. 

We are interested to see whether any diagnostic 
evidence exists for the greater predictability of 
blocking found by Bengtsson. Although our sample 
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sizes are fairly small, several interesting features 
related to differences in the nonlinear interaction 
terms appear. 


2. Procedures 

2.1. Data 

The data used in the study are the output of the 
twice daily operational analysis of the National 
Meteorological Center (NMC) from the winters 
(DJF) of 1976-77 and 1978-79. Data for the 
horizontal wind (u, v\ geopotential height (z) and 
temperature (T) on a 2.5° x 2.5° grid at the ten 
mandatory levels in the troposphere (1000, 850, 
700, 500, 400, 300, 250, 200, 150 and 100 mb) 
were used. Vertical velocities calculated from the 
quasi-geostrophic w equation (assuming w = 0 at 
1000 and 100 mb) were used in the energetics 
calculation (Chen et o/., 1981). Fourier coefficients 
of a, t\ z, T and to were computed for every 2.5° 
of latitude at the 10 levels with the wavenumber 
expansion of these variables truncated after zonal 
wavenumber 18. 

During 1976-77, NMC employed the analysis 
scheme of Flattery (1971) in which only the 
rotational component of the wind is retained in 
u and v . In 1978-79, a multivariate, optimum 
interpolation analysis scheme was employed (Berg- 
man, 1979). In this case a small divergent com- 
ponent of the horizontal wind is retained. How- 
ever, the kinetic energy content of this divergent 
wind is quite small (Chen and Tribbia, 1981) and 
we do not feel that it will adversely affect the 
evaluation of the vorticity advection terms in the 
enstrophy budget calculation (see below). 

2.2 Analysts 

In the interest of brevity, the spectral energetics 
and spectral enstrophy equations will be presented 
in symbolic form only. The formulation of the 
spectra] energetics equations can be found in 
Saltzman (1957, 1970) and the formulation of the 
spectral enstrophy equations can be found in 
Chen and Tribbia (1980). Note that derivation 
of the enstrophy equations requires that we assume 
the wind field to be nondivergent. The 1978-79 
NMC wind field is not nondivergent, as mentioned 
earlier, but we feel the advection of vorticity 
by the divergent component of the wind will be 
small compared to advection by the rotational 


component and, therefore, the inclusion of a small 
divergent component in the horizontal wind will 
not affect our calculation. 

The rates of change of wavenumber m kinetic 
energy (K m ), available potential energy (A m ) and 
enstrophy (E m ) are given by: 

S 

— C(tf „ KJ + C(A m . + C k (m\n, /) 

+*(*.) 

— A m — C(A iy A m ) ~~ C(A m , K m ) + C 4 (mtn, /) 

+ GWJ+JWJ 

and 

E m = C(£„ EJ + C f (mln, /) + fi m + G(E m ) 

+ B(E m ). 

Here, a z subscript denotes a zonal mean quantity 

l() 2 = — f ()<U where X = longitude], 

In Jo 

and an m subscript denotes a departure from this 
mean for wavenumber m. The notation C(A y B) 
represents a conversion of energy (or enstrophy) 
from reservoir A to reservoir B. C x (min, /) denotes 
the rate of increase of the quantity X (where X = K y 
A or £) at wavenumber m due to nonlinear triad 
interactions with all possible combinations of wave- 
numbers n and /. fi m represents the so-called 
^-effect (—V • Vf) in the enstrophy budget. We 
were unable to calculate the generation terms 
(GOO and G(E m )) with available data. B(X) 
denotes boundary flux terms due to fluxes of X 
across both vertical and horizontal boundaries. 

We integrated the budget equations over the 
region from 30° N to 80° N and from 1000 mb to 
100 mb. When not integrated over the entire mass 
of atmosphere, the budget equations include certain 
horizontal and vertical flux terms (Saltzman, 1970; 
Chen and Tribbia, 1980). With our data, we were 
unable to calculate the vertical flux terms, al- 
though the vertical fluxes could in principle, be 
competitive with the other terms. The horizontal 
flux terms are frequently small (except the flux of 
geopotential energy) and differences in them 
between blocking and nonblocking situations were 
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usually not significant. Therefore, we will not 
discuss them. 


2.3. Selection of blocking days 
Traditionally, blocking has been identified sub- 
jectively according to certain kinematic properties 
of the flow as set down by Rex (1950). Recently, 
blocking events have been identified as persistent 
departures from the climatological mean height 
field (Dole, 1978; Chamey et al., 1981). Our 
primary goal is to study persistent, large-scale 
blocking events compared to predominantly zonal 
circulations and to eliminate short duration (syn- 
optic time scale), low amplitude, small-scale 
features from our sample of blocking observations. 
Therefore, we restrict ourselves to cases of station- 
ary or slowly propagating ridges where the de- 
parture of the 500 mb height from the zonal mean, 
averaged over every 2.5° of latitude from 55°- 
80° N, exceeded 250 m for seven days or more. 
Blocking events so selected were compared with the 
twice-daily NMC 500 mb charts. The corre- 
spondence between the dates of blocking events 
identified subjectively from the charts and those 
determined from persistent height departures was 
quite good. Discrepancies of one or two ob- 
servations (12-24 hours) at most occurred at the 
beginning and ending of the individual blocking 
events. Frequently, the two selections were identi- 
cal. Therefore the blocking days identified from the 
height departures were used to form the sample of 
blocking and nonblocking days. Any observation 
not falling in a blocking period was included in the 
nonblocking sample. We should note that chang- 


ing one or two observations at the beginning or 
ending of a given blocking event from the blocking 
sample to the nonblocking sample (or vice versa) 
did not significantly affect the mean energy and 
enstrophy budgets. 

In addition, a case of very large negative height 
departures occurred during January and February 
1977. This particular event has been discussed in a 
number of contexts by numerous other authors. Its 
existence is correlated with a simultaneous, very 
cold sea surface temperature anomaly in the North 
Pacific Ocean (Miyakoda and Rosati, 1982). 
Chamey et al. (1981) provide the dates of the two 
periods during January and February 1977 of the 
largest negative height departures from the climato- 
logical mean. We will use their dates for the 
persistent low height anomalies in our diagnostic 
calculation. In fact, these low anomalies persisted 
throughout most of January and February 1977 
(Chen and Shukla, 1983). Because of the existence 
of large amplitude features of one type or another 
throughout the 1976-77 winter, no nonblocking 
days from this winter were included in our 
nonblocking sample. A summary of the blocking 
and nonblocking days for the two winters is given 
in Table 1. 

Our selection procedure is admittedly arbitrary 
but it does identify periods of large persistent height 
departures. The comparison of these periods with 
synoptic charts indicates that these periods satisfy 
(for the most part) the conventional, subjective 
definition of blocking. However, the inherent 
subjectivity in identifying what does or does not 
constitute blocking continues to be a problem in 
diagnostic studies of the phenomenon. 


Table 1. Tabulation of blocking and nonblocking days from the 1978-79 and 1976-77 winters determined 
from persistent height departures ( see text). Dates for the negative height anomalies in January-February 
1977 were determined by Chamey et al. (1981) as extreme departures from the climatological mean 
height 


1978-79 

Blocking days (90 observations) 

1. 0000 GMT 1 Dec.- 1200 GMT 7 Dec.: Atlantic 

2. 0000 GMT 20 Dec.- 1200 GMT 27 Dec.: Atlantic 

3. 1200 GMT 29 Dec.-OOOO GMT 10 Jan.: Pacific 

4. 0000 GMT 14 Jan.-1200 GMT 26 Jan.: Atlantic 

5. 0000 GMT 16 Feb.-0000 GMT 22 Feb.: Atlantic 
1976-77 

Blocking days (36 observations) 

1. 1200 GMT 1 1 Dec.-OOOO GMT 29 Dec.: Atlantic 


Nonblocking days (85 observations) 

0000 GMT 8 Dec.- 1200 GMT 19 Dec. 

0000 GMT 28 Dec.- 1200 GMT 28 Dec. 

1200 GMT 10 Jan.-1200 GMT 13 Jan. 

0000 GMT 27 Jan.- 1200 GMT 15 Feb. 

1200 GMT 22 Feb.-1200 GMT 28 Feb. 

Chamey et al. (1981) negative anomalies (76 observations) 

1. 30 December-1 7 January 

2. 2 February-22 February 


TeUus 36 A (1984), 1 





ORIGINAL PAGE IS 
QE POOR QUALITY 


A COMPARISON OF THE SPECTRAL ENERGY AND ENSTROPHY BUDGETS 55 


3. Results 

3.1. 1978-79 winter 

From December 1978 through February 1979, 
we identified five major blocking events, four over 
the Atlantic-European region and one over the 
North Pacific (Table in general, the mean 
energy and enstrophy budgets of the Atlantic cases 
were similar to the Pacific case even though the 
initiation of the Pacific case was different from the 
late December 1978 case (Hansen and Chen, 
1982), Whether or not the energetics of the Pacific 
case are typical of other Pacific blocking cases is 
unknown, although all the Atlantic events showed 
similar characteristics and were similar to the 
December 1976 case (see below). First, consider 
the average kinetic energy {K\ enstrophy (£) and 
available potential energy (A) spectra (Fig. 1) of all 
five cases compared to the nonblocking days. In the 
mean, the kinetic energy in the longest waves 




Fig. 1 . Wavenumber spectra of (a) the kinetic energy, 
K m \ (b) the enstrophy, E m \ and (c) the available potential 
energy, A mf for the 1978-79 winter. The solid line 
represents the mean spectra for all blocking days and the 
dashed line represents the nonblocking days. 


C CA m . K m > 





Fig . 2. 1978-79 energy budget terms: (a) CM If A J; (b) C{A m , it*); (c) C(K t > K m )\ and (d) C*(mln, /)• As in Fig. 1, 
the solid line signifies the spectra for all blocking days and the dashed line denotes all nonblocking days. The notation 
is explained in the text 
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(wavenumbers 2-3) is greater during blocking 
compared to the nonblocking sample, while the K 
content of the intermediate-scale wavenumbers is 
slightly lower during blocking. The enstrophy 
sp e ctru m is qualitatively similar to the K spectrum 
at low wavenumbers but the total enstrophy in 
wavenumbers 5 through 10 is nearly the same for 
blocking and nonblocking. The A spectrum shows 
only an increase in available potential energy at 



Fig. 3. C k (m\n % f) for the four Atlantic cases from 
1978-79 (solid line) compared to all nonblocking days 
(dashed line). 


wavenumber 1 for blocking and no other sys- 
tematic differences. 

The energy conversions C(A X% A m ) and 
C(A m , K m ) (Fig. 2a, b) are quite similar for the 
blocking and nonblocking samples. C(A s ,A m ) 
shows slight increases at wavenumbers 1 and 3 
during blocking with largely offsetting differences 
in wavenumbers 4 through 7. The increased 
conversion at wavenumbers 1 and 3 is totally 
absent when the four Atlantic cases are considered 
alone (not shown). The C(A m ,K m ) spectrum is 
remarkably similar in the blocking and nonblock- 
ing situations although there are slightly larger 
conversions in wavenumbers 7 through 10 during 
blocking. C(K l% K m ) (Fig. 2c) indicates a large gain 
in wavenumber 3 kinetic energy due to mean flow 
interaction that is largely offset by reduced values 
at wavenumbers 1 and 2 when comparing blocking 
and nonblocking episodes. 

The largest systematic differences between 
blocking and nonblocking appears in the nonlinear 
interaction term, C*(mln, /) (Fig. 2d). The transfer 
of kinetic energy from intermediate-scale wave- 
numbers to the largest-scale wavenumbers takes 
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Fig. 4. The 1978-79 enstrophy budget terms: (a) C^min, /); (b) C(£ x , E m ); and (c) f m% the rate of change of E m due 
to the 0 effect As in Fig. 1, the solid line represents all blocking days in the 1978-79 winter and the nonblocking 
results are represented by the dashed line. 


Tellus 36 A (1984), 


A COMPARISON OF THE SPECTRAL ENERGY AND ENSTROPHY BUDGETS 


57 


place at a much greater rate during blocking 
periods compared to nonblocking, mostly into 
wavenumber 1. For the Atlantic cases alone (Fig. 

3) , both wavenumbers 1 and 3 benefit from the 
much increased kinetic energy export from wave- 
numbers 6 through 9 during blocking. C A (m\n, f) 
(not shown) is similar for the two samples and the 
most important boundary flux term was the 
horizontal flux of geopotential energy (not shown) 
which contributed more kinetic energy to wave- 
number 1 during blocking with an offsetting loss at 
wavenumber 3. 

The dominant term in the enstrophy budget (Fig. 

4) is the nonlinear exchange term, C f (ml«, /)• A 
very remarkable difference between blocking and 
nonblocking situations appears. Normally, ens- 
trophy is cascaded from the longest waves into the 
shortest wavelengths (Steinberg et al., 1971). This 
characteristic also appears in our nonblocking 
sample (Fig. 4a, dashed line). However, during 
blocking, the net enstrophy cascade out of the 
longest waves (wavenumbers 1-5) disappears and 
a net gain of enstrophy appears at wavenumber 2. 
Meanwhile, the intermediate wavelengths (wave- 
numbers 6-10) are exporting a significantly larger 
amount of enstrophy than during nonblocked 
periods. For the Atlantic cases alone (not shown) 
there is an increase in enstrophy due to nonlinear 
interactions at both wavenumbers 2 and 3 during 
blocking. Thus during blocking the enstrophy 
cascade at low wavenumbers is reversed , as 
enstrophy is injected into the largest-scale waves 
by transient cyclone-scale waves. The only other 
difference in the enstrophy budgets is a greater 


loss of wavenumber 3 enstrophy due to the effect 
during blocking. 

The differences in the nonlinear interaction terms 
are easily seen by comparing the kinetic energy and 
enstrophy flux functions for our two samples. 
Following Steinberg et al. (1971), these functions 
can be defined as 


El 

6m 


— C k (m\n,f) 


and 


El 

6m 


— C £ (mln, l) 


assuming £*(0) = F B ( 0) = 0. Regions of the ap- 
propriate flux function with negative slopes re- 
present wavenumber bands gaining AT or £ through 
nonlinear interactions. Regions with positive slopes 
are losing energy or enstrophy. The absence of an 
enstrophy flux out of the lowest wavenumbers 
during blocking and the change in the sign of the 
enstrophy flux function near wavenumber 2 is 
evident in Fig. 5a. The much larger flux of kinetic 
energy into the long waves during blocking can 
clearly be seen in Fig. 5b. The flux functions for the 
Atlantic cases alone are essentially the same as 
those shown in Fig. 5. 


3.2. 1976-77 winter 

A comparison between the energy and ens- 
trophy budgets of the December 1976 case of 
Atlantic blocking and the persistent, large am- 
plitude low anomalies during January and Feb- 




Fig. 5. The nonlinear flux functions and (a) enstrophy, £ e , and (b) kinetic energy, £*, for blocking (solid line) and 
nonblocking (dashed line) for the 1978-79 winter. 
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ruary 1977 arc given in Figs. 6-8. In general, the 
December 1976 event showed relatively weak 
large-scale flow wavenumber) baroclinic processes 
compared to January-February 1977 (Fig. 7a, b) 
(or the 1978-79 nonblocking sample), but it does 
exhibit the same enhanced upscale kinetic energy 
flux as the 1978-79 blocking cases (Fig. 7d), and a 
very striking reversal of the low wavenumber 
enstrophy cascade (Fig. 8a) at wavenumber 2. The 
similarity in the behavior of Rex -type blocking in 
two different winters lends support to the sig- 
nificance of this result. 

However, comparing the December 1976 event 
to the January and February 1977 events indicates 
remarkable changes in the circulation between 
December 1976 and the following two months. 
There was a remarkable increase in C(A Z% A m ) and 



C(A m , K m ) at wavenumbers 2 and 3 during 
January-February 1977 (Fig. 7a, b). This char- 
acteristic was also noted by Chen and Shukla 
(1983) but is even more striking when compared to 
the December 1976 statistics. There is also a larger 
net gain at low wavenumbers in K m due to 
C(K X , K m ) in January-February 1977 (Fig. 7c). 
The impression given by these results is that 
planetary-scale baroclinic and barotropic pro- 
cesses were of dominant importance in the 
January-February 1977 case as opposed to the 
case for Rex-type blocking. The nonlinear, wave- 
wave interactions during January-February 1977 
were dominated by large energy gains for wave- 
numbers 1, 3 and 4 and a large loss from 
wavenumber 2. The kinetic energy transfer out of 
the intermediate-scale waves was relatively small. 




Fig. 6. Wavenumber spectra of (a) K m \ (b) E m \ and (c) A m for the 1976-77 winter. Here the solid line signifies the 
December 1976 Atlantic blocking case and the dashed line represents the periods of large negative height anomalies in 
January and February 1977. 
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Fig. 7. (a) C{A l% A m )\ (b) C(A m . tf w ,): (c) CC*,. KJi and (d) C*(mtn, /) for the 1976-77 winter. As in Fig. 6, the 
solid line represents the December 1976 blocking case and the January-February 1977 negative anomalies energetics 
are represented by the dashed line. 


This is consistent with the reduced intermedi- 
ate-scale baroclinic activity | C(A t ,A m ) and 
C(A m , K m ) at wavenumbers 7-10 and 6-10, 
respectively (Fig. 7a, b)l during January-February 
1977. 

The nonlinear enstrophy exchange, C £ (mln, /) 
(Fig. 8a), for the December 1976 case is similar to 
the 1978-79 cases as noted above (e.g., reversal of 
cascade at low wavenumbers and increased export 
of enstrophy from intermediate-scale wave- 
numbers), but the January-February 1977 ens- 
trophy exchange is more like the 1978-79 non- 
blocking sample. These qualitative features carry 
over to the enstrophy and energy flux functions 
(Fig. 9). Notice the very striking reversal in the 
enstrophy flux function at low wavenumbers for the 
December 1976 case (Fig. 9a). Also notice the 
much lower value of the enstrophy flux function for 

Tellus 36A (1984). 1 


the January-February 1977 case in a nearly 
constant range from roughly wavenumbers 8 to 15. 
The enstrophy flux function in the same range for 
the December 1976 blocking case is virtually 
identical to the 1978-79 blocking and nonblocking 
samples. 

3.3. Predictability question 
The rate of error growth in numerical weather 
prediction models due to small-scale inaccuracies in 
the initialization can be related to the rate of 
down-scale enstrophy cascade in a constant ens- 
trophy flux inertial subrange of a two-dimensional 
turbulent fluid (Leith and Kraichman, 1972). The 
characteristic time of this process can be called the 
predictability time and is inversely proportional to 
the cube root of the constant rate of enstrophy 
transfer to higher wavenumbers in a constant 
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Fig. 8. The 1976-77 enstrophy budget terms: (a) C f (mln, /); (b) C(£,» £*); and (c) 0 m for the December 1976 
blocking case (solid line) and the January-Febniary 1977 negative anomaly cases (dashed line). 



wavanumbar aovanumbar 

Fig. 9 . The nonlinear flux functions for (a) enstrophy and (b) kinetic energy for the December 1976 blocking case 
(solid line) and the January-Febniary 1977 negative anomalies (dashed line). 


enstrophy flux inertial subrange (Lflly, 1970; Leith bert, 1981), the fact that the atmosphere at the 

and Kraichnan, 1972). Although efforts to find larger scales behaves predominantly like a two- 

observational evidence of a two-dimensional tur- dimensional fluid in which the kinetic energy 

bulent subrange in the atmosphere have been spectrum obeys a —3 power law (Chen and 

inconclusive (Chen and Wiin-Nidsen, 1978; Lam- Tribbia, 1981; Julian et al., 1970; Wiin-Nielsen, 
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1967) suggests that a qualitative estimate of the 
predictability time for blocking compared to 
nonblocking situations based on the enstrophy flux 
function may be made. The estimate can only be 
done in a qualitative way because the numerical 
value of the enstrophy flux function, F B% in the 
constant range is dependent upon the truncation 
used (Steinberg et al., 1971). The reason for this 
sensitivity is readily apparent from inspection of the 
enstrophy spectra (Fig. lb or Fig. 6b). A finite 
amount of enstrophy is lost when the wavenumber 
expansion of the wind field is truncated at wave- 
number 18. This will cause aliasing problems in the 
nonlinear enstrophy interaction in our truncated 
data set that will cause changes in the shape of the 
enstrophy flux function at high wavenumbers. 
Unfortunately, the accuracy of our data at wave- 
numbers higher than 18 is poor so this problem 
cannot be alleviated. However, we will assume the 
values of the flux function that we have are 
qualitatively correct for comparison purposes. 

As noted earlier, the enstrophy flux function in 
the constant range is nearly identical in both 
blocking and nonblocking situations (Figs. 5a and 
9a] suggesting no greater predictability for block- 
ing. However, the January-February 1977 event, 
which was characterized by enhanced large-scale 
baroclinic processes, exhibited a roughly 1/3 
smaller enstrophy flux rate (Fig. 9a, dashed line). 
This result allows us to speculate that the 
persistent low anomalies in 1977 were inherently 
more predictable than either the nonblocking or 
Rex-blocking cases. 

An explanation of the greater predictability of 
the Rex-type blocking found by Bengtsson (1981) 
might be found in the reversal of low wavenumber 
enstrophy cascade (especially near wavenumber 2), 
and the absence of an enstrophy cascade out of the 
lowest 5 wavenumbers during blocking (Figs. 5a 
and 8a) as opposed to the normal case (Fig. 5a 
(dashed line) and Steinberg et al., 1971). This may 
indicate that the destruction of large-scale vorticity 
by down-scale cascade and eventual dissipation 
is greatly reduced or eliminated during blocking 
events. As a result, blocking patterns may be more 
persistent and therefore more predictable. 

We should re-emphasize that the source of the 
enstrophy being cascaded into the low wave- 
numbers is the intermediate-scale wavenumbers. 
The fact that the predictability time for blocking is 
the same as for nonblocking may indicate that 


blocking is no more predictable in the traditional 
sense than the cyclones that successively re- 
generate it. 

K 

4. Concluding remarks 

To summarize, the most noteworthy differences 
between the spectral energy and enstrophy budgets 
of blocking versus nonblocking periods are in the 
nonlinear interaction terms. Apparently, when the 
ultralong waves have a large amplitude, such as 
during Rex-type blocking, they are able to more 
efficiently gain kinetic energy from the cyclone 
scale waves through wave-wave interactions. In 
addition, through a reversal of the low wave- 
number enstrophy cascade, the ultralong waves are 
also able to gain enstrophy from the intermedi- 
ate-scale waves during blocking, whereas the 
ultralong waves lose enstrophy through wave-wave 
interactions in nonblocked circumstances. The 
absence of enhanced large-scale baroclinic pro- 
cesses during the 1978-79 and December 1976 
Rex-type blocking events does not augur well for 
the applicability to these cases of the various 
baroclinic blocking theories proposed recently. 

The fact that the energy and enstrophy budgets 
for Rex-type blocking were similar in the two 
winters studied suggests that our findings are 
typical. However, our results for changes in the 
nonlinear interaction terms should be tested with a 
larger data set. This calculation is currently in 
progress. 

The January-February 1977 negative height 
anomaly case presents a different class of pheno- 
menon in which low wavenumber baroclinic energy 
conversions maintain the large amplitude features 
of the circulation and in which a significantly lower 
enstrophy flux function suggests greater predict- 
ability. This is in direct contrast to the Rex-type 
blocking cases which show no greater predict- 
ability in the traditional sense compared to non- 
locking. However, the reversal of the low wave- 
number enstrophy cascade does suggest a re- 
duction in the dissipation of the low wavenumber 
enstrophy (assuming the atmosphere behaves like a 
two-dimensional turbulent fluid) and therefore 
greater persistence (and predictability) of the 
blocking pattern. In addition, the enstrophy cas- 
cade reversal may be an indication of intermittent 
behavior (Batchelor, 1960) in the large-scale flow 
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requiring no particular low wavenumber in- 
stability. (This aspect will be discussed further in a 
forthcoming paper.) If so, a description of blocking 
in terms of vortex dynamics (Batchelor, 1960; 
Onsager, 1949) may be useful. This possibility is 
being studied further. 

From our results it appears that cyclone-scale 
waves play an important role in maintaining 
quasistationary blocking patterns by injecting 
enstrophy and kinetic energy into the large-scale 
flow. Bengtsson (1981) has also noted the ap- 
parent importance of transient smaller-scale waves 
in one case of blocking during January, 1979. 
Saltzman (1959) first suggested that nonlinear 
kinetic energy transfer from transient cyclone-scale 
eddies acts to maintain the quasistationary centers 
of action. Holopainen and Oort (1981) have 
suggested that in extratropical latitudes transient 
eddies act to maintain the enstrophy of the standing 
waves although more recent work disputes the 


earlier result (Holopainen et al., 1982). The role of 
transient eddies in maintaining quasistationary 
features of the circulation deserves further study. 
Also, it would appear that theoretical models of 
blocking should incorporate the effects of transient 
eddies. 
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ABSTRACT 

Average statistics for periods oflarge positive 300 mb height anomalies are compared to statistics for all other 
situations using NMC data for the 13 Januaries from 1963 to 1977. The 500 mb heights and geostrophic 
streamfunctions are represented as surface spherical harmonics, and energy and enstrophy spectra along with 
nonlinear wave-wave interaction statistics are computed. 

Differences in 300 mb geopotential height variance, kinetic energy and enstrophy spectra occur between large 
positive anomaly events and other days in the twodimensional spectral index band from roughly « - 6 to 
»*9, where n is the degree of the associated Legendre function. The same index band experiences a reversal 
of both the usual kinetic energy and enstrophy cascades during large positive anomaly events. That is, the 
6 < n < 9 band gains energy and enstrophy fiom wave-wave interactions during the anomaly events and loses 
energy and enstrophy by the same process at other times. The source of this energy aod enstrophy is higher 
index (smaller twodimensional scale) waves. The indication is that the Atlantic cases are more subject to this 
cascade reversal than are Pacific events. 

Our results suggest that the smaller scale, transient eddies may play a regime-dependent role in interactions 
with atmospheric circulation modes on the scale of the persistent anomalies. When interacting with larger-scale 
features, the role of smaller-scale transients may not always be dissipative. 


1. Introduction 

Studies of the interaction between transient eddies 
and the climatological-mean stationary flow indicate 
that although the enstrophy budget of stationary waves 
in midlatitudes appears to be maintained in part by 
the transient eddies (Holopainen and Oort, 1981), the 
effect of transient eddy heat transports overwhelms the 
vorticity transports in a potential vorticity budget, giv- 
ing an overall dissipative role to the transients when 
interacting with the time mean flow (Youngblut and 
Sasamori, 1980; Holopainen etal., 1982). Further study 
has revealed that the long period transients (periods of 
10 to 90 days) make a stronger contribution to this 
dissipative effect than do synoptic-scale (2.5 to 6 day 
period) transients (Lau and Holopainen, 1984). Van 
den Dool (1983), on the other hand, infers a potentially 
large role for high frequency transients in forcing 
monthly mean anomalies. 

Theoretical studies have indicated that the large-scale 
wave pattern can act to modulate synoptic-scale baro- 
clinic activity (Fredrickson, 1979) and that synoptic- 
scale transients can play an important role in main- 
taining equilibrium states, or “weather regimes'’, in 
certain models (Reinhold and Pierrehumbert, 1982; 
Kallen, 1981, 1982). In the present study, we will use 
observational data to examine whether the interaction 
between different scales of motion in a spherical har- 


monic framework is dependent upon the circulation 
regime characterizing the large-scale pattern. We use 
the large positive height anomaly events identified by 
Chamey et al. ( 1 98 1 ) to stratify the data in the 1 5 Jan- 
uaries from 1963 to 1977. Statistically significant dif- 
ferences in the energy and enstrophy spectra as well as 
in the nonlinear wave-wave interactions are found. 

Baer (1972) suggested that the degree of the Legendre 
function in a spherical harmonic representation of the 
data is the appropriate index to use when studying at- 
mospheric spectra. Observational studies of nonlinear 
exchanges using this two-dimensional index have been 
made by Chen and Wiin-Nielsen (1978), Lambert 
(1981), and Boer and Shepherd ( 1983). Mean cascading 
during winter as calculated by Chen and Wiin-Nielsen, 
for example, shows that for n • 1 to 7 and n = 19 to 
31, kinetic energy is gained through nonlinear inter- 
action, with this energy being supplied by the modes 
with n * 8 to 18. In the enstrophy budget, modes with 
n » 1 to 7 gain a very small amount of enstrophy while 
a large loss of enstrophy from the range n * 8 to 18 is 
balanced by gains of higher index modes. The physical 
explanation of the occurrence of these wave-wave in- 
teractions in the atmosphere is that the dissipation of 
kinetic energy and enstrophy takes place in different 
wavenumber bands from where these quantities are 
generated. In the present study, a dataset spanning a 
longer time period than those considered in earlier 


e 1986 American Meteorological Society 


I March 1986 


BROWN, JR., HANSEN AND PANDOLFO 


477 


studies is used and attention is focused on interactions 
near the low wavenumber end of the spectrum. 

First, we illustrate the differences in the height vari- 
ance, kinetic energy, and enstrophy spectra of the large 
positive anomaly events (LPAs) compared to more 
normal circulations. The changes occur primarily in 
the spectra] index band from n = 6 to 9. Guided by 
these results, we explore systematic differences in the 
nonlinear cascading. We find an apparently regime- 
dependent role of the smaller scale transient eddies 
when interacting with the spectral band that distin- 
guishes the LPAs. 

2. Data 

The basic data used in this study are 500 mb geo- 
potential heights from the NMC final analyses. For 
each January day in the years 1963-77, daily values 
of the geostrophic stream function are computed. By 
using January data only, we hope to avoid any com- 
plications due to the seasonal cycle. We present results 
for the 500 mb surface only because these results may 
be considered representative of the vertical average in 
the energy and enstrophy budgets. (Qualitatively sim- 
ilar results were obtained using 300 mb data.) 

To perform the required spherical harmonic anal- 
yses, data are required over the entire sphere. The NMC 
grid, however, extends only from 90°N to about 1 5°N. 
Data must be extrapolated to fill the region from 15°N 
to the equator about which a symmetry assumption 
can be made to account for the Southern Hemisphere. 
In order to generate the required data, Ellsaesser's 
(1966) technique has been adopted. To construct values 
in the north equatorial region, data lying within the 
NMC octagon are used to extrapolate southward along 
the meridians using a linear function of sin(2<$) as an 
extrapolation profile (<t> *= latitude). Even symmetry of 
the geopotential height about the equator is then as- 
sumed. 

The fictitious data, of course, affect the spectral coef- 
ficients. Nevertheless, we would not expect the intro- 
duction of such data to keep us from our main objective 
of identifying essential differences between spectra for 
different periods as long as the manufactured data are 
based on reasonable, though crude, physical assump- 
tions. Comparison of our results obtained in this way 
to those of Boer and Shepherd (1983) in which global 
data were used reveals quite similar results. Therefore, 
we do not expect the treatment of the tropics in the 
present study to dramatically affect the results. 

In order to define periods of large positive anomaly 
events, we have adopted the criterion of Chantey et al. 
(1981). In their work, time-longitude cross sections of 
departures of 500 mb heights from their climatological 
mean values at 50 e N were examined for the 15 winter 
seasons that include the January months selected for 
this study. A period was designated as a large positive 
anomaly event if at any point on the 50°N latitude 


circle a positive anomaly of 200 gpm or more persists 
for 7 days or more. Some subjectivity was introduced 
in determining the beginning and end of each event 
(Charney et al., 1 98 1 ). Table 1 lists the periods consid- 
ered to be LPAs. The letter A, P or B in this table 
indicates whether the largest amplitude ridge occurred 
over the Atlantic or Europe (A), the Pacific Basin (P), 
or whether comparable amplitude ridges occurred over 
both oceans (B). Events listed of less than 7 day du- 
ration spanned the beginning or ending of the month 
in question. Note that the existence of a large positive 
height anomaly doesn’t necessarily imply that a greater 
height variance overall will exist at any particular lo- 
cation or for any particular wavenumber (e.g., see the 
discussion in Charney et al., 1981). 

Of the 465 days included in our analysis, there were 
154 LPA days and 31 1 “normal” days according to 
the criterion used. A comparison of the characteristics 
of blocking and nonblocking days determined from 
purely subjective review of individual twice-daily 
weather charts for the 15 months in question yielded 
generally similar results to those from the objectively 
chosen LPA cases. However, our purpose is to compare 
diagnostic statistics from time periods characterized by 
persistent circulation anomalies to statistics from time 
periods in which such features were lacking. The re- 
lationship of LPAs to blocking is outside the scope of 
the present study. 

3. Analysis 

We compute a streamfunction f from the geo- 
strophic balance equation: 

/Vty + V/-V*«gV 2 z (1) 

where /is the Coriolis parameter, g the acceleration of 
gravity, and z the height of the constant pressure sur- 
face. For the nondiveigent component of the flow, the 
kinetic energy ( K ) and enstrophy (£) can then be writ- 
ten as 


Table 1. January days (1963-77) on which 
LPA events took place. 


1 

8-21 January 1963 (P) 

2 

22-29 January 1963 (B) 

3 

30-31 January 1963 (P) 

4 

1-6 January 1965 (P) 

5 

1 January 1966 (P) 

6 

5-14 January 1967 (A) 

7 

1-5 January 1968 (B) 

8 

6-31 January 1968 (A) 

9 

1-2 January 1969 (B) 

10 

3-14 January 1969 (P) 

11 

1 1-25 January 1970 (A) 

12 

7-13 January 1971 (B) 

13 

15-3! January 1971 (P) 

14 

3-9 January 1973 (A) 

15 

1-10 January 1974 (P) 

16 

27-31 January 1975 (P) 

17 

6-12 January 1977 (A) 


478 


JOURNAL OF THE ATMOSPHERIC SCIENCES 


Vou 43, NO. 5 


ij.,1 



(2) 


(3) 


where m = sin<£ and 4> is latitude. Other notation is 
conventional. 

On a sphere, any real-valued function such as ^ can 
be written formally as a series of surface spherical har- 
monics in the form 


iKX, m) = 2 2 (A m m cos mX + B m m snrnX)P m m (ji) 

Jr— m 

(4) 

where P n m (n) is the associated Legendre function. The 
series is uniformly convergent (hence term-by-term dif- 
ferentiable) provided that it has continuous second de- 
rivatives over the whole sphere (Hobson, 1931, p. 602). 
Using the orthogonality properties of the Legendre 
functions, the coefficients A” and B” can be obtained 
from 



4AX, M){ C0Sm *}/V"(M)dM4 
IsinmX J 


(5) 


where S„ - 1 for m * 0, and S m = 2 for m + 0. 

Series differentiability is an important consideration 
since we are to examine spectra obtained by application 
of various differential operators to series representations 
of the streamfunction. In particular, since f = Vfy, 
differentiability of (4) allows us to express the vorticity 
as 


flX, m) - 2 2 ~ ^ (A„ m cos mX 

m-0 i*- iff u 

+ B m m sinmXJ/W) (6) 

where a is the earth's radius. 

The kinetic energy and enstrophy can then be written 
in series form as 


K = 2 K„ = 5 2 n( ”^ — } i m [{Arf + (BiTf] 

it -0 n -0 m -0 ^ 

(7) 

cd ® * - 2 /- 4 * 

£ - 2 - 2 2 Wr? + Wfl. 

w -0 »-0 m -0 ^ 

( 8 ) 

Equations governing the complete spectral energy 
and enstrophy budgets of the atmosphere can be ob- 
tained by differentiating (2) and (3) with respect to time 
and substituting from the vorticity equation. The ap- 
propriate vorticity equation for the atmosphere in- 


cludes generation and dissipation terms as well as terms 
representing the advection of relative and planetary 
vorticity. It isn’t possible to reliably estimate the former 
two terms with the available data so we will focus at- 
tention on the barotropic terms we can compute. 
Therefore, we won’t present complete budget statistics 
but will emphasize a particular aspect of the problem, 
namely barotropic, nonlinear effects. The nondiver- 
gent, adiabatic, frictionless form of the vorticity equa- 
tion may be written as 


dt 


-JW, D - 


2ft d* 

a 2 ax’ 


(9) 


where J denotes the Jacobian 


f) 


_d^d£\ 

a 2 W dX dX dj' 


( 10 ) 


We can obtain simplified spectral forms of the kinetic 
energy and enstrophy budget equations by differen- 
tiating (2) and (3) with respect to time and using (4), 
(6) and (9). In particular, the derivatives of K and E 
can be represented as functions of n as: 



n 


2 i» m 


= 2 jwW + W. (li) 

n 4 


2 J. 

m-0 


= 2 \ ° + B”rr>. (12) 

m-0 ^ a 

Here, q„ m and r„ m are the amplitudes of the cosine and 
sine components of the harmonic expansion of J(\fr, 
f). That is. 



M, (13) 

IsinmX ) 


I” and J” measure the rate of gain (loss) of kinetic 
energy and enstrophy respectively in component (n, 
m) that results from a corresponding rate of loss (gain) 
in the remainder of the spectrum. The change in total 
kinetic energy or enstrophy due to wave-wave inter- 
actions is zero. Computational considerations related 
to calculating these spectra and nonlinear transfers are 
discussed in the Appendix. 

As noted earlier, for a less approximate form of the 
vorticity equation, (II) and (12) would also contain 
generation and dissipation terms. When averaging over 
a large ensemble of days, the tendency term in the en- 
ergy or enstrophy budget will effectively be zero so the 
effects of the wave-wave interactions must be balanced 
by the effects of generation and dissipation for a given 
value of n. 

In truncating (4) or (6), we wish to retain all terms 
that correspond to waves resolvable by the original data. 
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ir it is assumed that four grid points are required to 
resolve a wave adequately, the limit of resolution of 
our data is near zonal wavenumber m = 18. With 
comparable data reliability in both latitude and lon- 
gitude, the limit of resolution in latitude as measured 
by the one-dimensional index n — m (the number of 
zeroes of P”(p) between the poles) is also 1 8. However, 
for the sake of easy comparison with earlier studies, 
we have employed triangular truncation at n = m 

- 30 in order to retain information on higher order 
analysis modes that are present (e.g., Chen and Wiin- 
Nielsen, 1978). The effect of differing truncations on 
the nonlinear interaction statistics is also discussed in 
the Appendix. The results converge for increasing 
truncation even though the values of or J„ m are, in 
general, a function of the truncation (see Appendix). 
The values of /„ or J n appear to be insensitive to the 
truncation for n *£ 14. 

4 . Results 

a. Spectra 

Austin (1980) has shown that blocking situations 
are characterized by large amplitudes and normal 
phases for the zonal harmonic wavenumbers m = 1, 
2, and 3. Pacific blocks in the mean exhibit greater 
amplitude for m - 2 and 3 while Atlantic blocks nor- 
mally exhibit greater amplitude in wavenumbers 1 
and 2. 

To illustrate the characteristic signature of LPAs in 
the two-dimensional index n, let us first consider the 
spectrum of geopotential height variance, 

(Z zr); ( 14 ) 

Z„ m , the amplitude of spherical harmonic (n, m ) for 
the 500 mb geopotential height, follows straightfor- 
wardly from the expansion of z in a series similar to 
(4). Figure 1 shows that the largest difference in a 2 
between LPAs and the balance of the data occurs at n 
= 7, 8 where LPAs have a greater variance. The LPAs 
also have slightly greater variance at n * 5 and slightly 
less at n * 3. 

Examination of K and E spectra brings out an index 
band where differences between LPAs and the remain- 
der of the sample are more systematic. The LPAs have 

greater kinetic energy in roughly the range 6 < n < 9 

• 

(Fig. 2). The average increase in K m is, 30 percent 

«-6 

for the LPAs compared to all other days. The statistical 
significance of this increase was tested with a two-sided 
/-test (Panofsky and Brier, 1958) assuming 5 days be- 
tween independent samples of the 500 mb height. The 
differences in the K H spectra are significant with greater 
than 99 percent confidence at r? = 7, 95 percent at n 

- 8, and 90 percent at n * 6. The differences in K„ for 
the Atlantic cases alone are greatest at n * 6, 7 while 



FJC. 1. Geopotential height variance a, 2 versus n 
for the 500 mb level. 


the Pacific cases show particularly prominent differ- 
ences at n ■ 7, 8 with smaller differences at n = 6 and 
9 (not shown). Differences in the enstrophy spectrum 
(Fig. 3) closely follow those in K n since E„ - [n(n 
+ 1 )/a 2 )K n . The stationary contribution to the K„ and 
£„ spectra was computed from the 15 January mean 
streamfunctions and is also displayed in Figs. 2 and 3. 
Note that the climatological mean stationary waves 
make significant contributions to the total K„, for ex- 
ample, only for n < 6. For n & 10, the spectrum is 
dominated by transient waves (see also Boer and Shep- 
herd, 1983). 

Larger values of K n or E n imply greater dissipation 
since the dissipation is most likely related to the inten- 
sity of the wind in some way. Therefore, we might ask 
what mechanisms balance this dissipation and allow 
the noted increases in K„ and whether smaller scale 
waves act in a dissipative role. 

b. Nonlinear exchange 

We must keep in mind the inherent artificiality of 
describing atmospheric fluctuations in terms of the 
chosen (or any other) basis functions. The results for 
individual wavenumbers are often difficult to interpret. 
Therefore, it is traditional to form ensembles of wave- 
numbers with which to synthesize desired components 
of the flow (Saltzman and Fleisher, 1960; Wiin-Nielsen 
etal., 1963; Blackmon, 1976; Hansen and Chen, 1982). 
By analogy with the basic properties of Fourier trans- 
forms, these ensembles of wavenumbers allow specific 
features of the circulation to be synthesized (e.g.„ Han- 
sen and Chen, 1982). 

Blackmon ( 1 976) has shown that geopotential height 
fluctuations with periods of 10-90 days are most im- 
portant in the wavenumber regime, 7 < n < 12. Lesser 
contributions to fluctuations of this period are made 
in the 1 < n < 6 range. This is consistent with our 
finding that the major difference in the kinetic energy 
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Fig. 2. Kinetic energy AT. versus n for the 500 mb level for LPA 
events (solid line); non- LPA events (dashed line) and for the 1 5 Jan- 
uary climatological mean stationary waves (dotted line). 


spectra between LPAs and the balance of our data occur 
in the range 6 < n < 9. 

The selection of the wavenumber ensembles we shall 
use presently was based on the differences in the spectra 
and was motivated as follows. First of all, harmonics 
with I < n < 5 are grouped together because these 
modes represent the gross, mean features of the cir- 
culation and are very steady in time (Blackmon, 1976; 
Boer and Shepherd, 1983). In the present dataset, the 
content of K„ and E„ for 1 < n < 5 is due almost 
entirely to the climatological mean waves (Fig. 2). In 
addition, there is virtually no difference in the mean 
K„ spectra, for example, between LPAs and the other 
situations for these modes (Fig. 2). 

A second ensemble is formed from the harmonics 
with 6 < n < 9. This grouping is due to the mean 
increase in the K and E content of this ensemble during 
LPA events. Note that this ensemble contains zonal 
wavenumbers from m - 0 through m = 9 and includes 
1 8 different spherical harmonics. Thus, depending on 
the relative amplitudes and phases of the components 
in specific realizations, relatively localized features can 
be described. A third ensemble corresponding to 
smaller two-dimensional scale waves is formed from 
waves with n > 10. These ensembles provide a con- 
venient way to visualize the results, but they are by no 
means unique. For example, based on differences in 
the cascading alone, ensembles for 1 < n < 3 and 4 
<n< 10 could be justified. 

The formation of wavenumber ensembles in the two- 
dimensional index is not directly analogous to the pro- 
cedure in the one-dimensional index used, for example, 
by Hansen and Chen ( 1 982). Each of our three chosen 
ensembles for the index n contains zonal harmonic 
wavenumbers of planetary scale (e.g., Blackmon, 1976). 
However, these waves have different latitudinal scales 
and therefore different two dimensional indices. 

The results of the calculation of /„ and J m are pre- 
sented in Tables 2 and 3, respectively. 


The gain in energy for the 1 < n < 3 ensemble is 
reduced during LPA events. This change is entirely 
due to the change at n = 3 which is partially offset by 
the changes at n = 4 and 3. Differences in the enstrophy 
cascade for this ensemble are small. However, we can 
see that in the range of 6 < n < 9, when all the positive 
anomaly cases are taken together, there is a reversal in 
the sign of the kinetic energy and enstrophy cascade 
compared to the balance of the sample. That is, during 
the positive anomaly events, waves with 6 < n < 9 gain 
both kinetic energy and enstrophy from the other waves 
due to wave-wave interactions; whereas at other times, 
this ensemble loses a comparable amount of K„ and 
E h due to the same mechanism. The source of this 
energy and enstrophy is the higher n waves, especially 
those for which 10 < n < 14, a region of the spectrum 
dominated by transient eddies (e.g.. Fig. 2 of Boer and 
Shepherd, 1983). The difference in the kinetic energy 
cascade between LPA and non-LPA days. A/*, is il- 
lustrated in Fig. 4. The implication is that smaller two- 
dimensional scale transients are making positive con- 
tributions to both the energy and enstrophy budgets of 
the LPA pattern during these events. This is a some- 
what remarkable result considering the global-scale 
treatment of the data used here. These results represent 
average interactions over the entire Northern Hemi- 
sphere. 

When the LPA events are subdivided into Atlantic 
and Pacific cases, the reversal of the cascades for the 6 
< n < 9 regime is much more pronounced for the 
Atlantic cases. For the Pacific cases, the enstrophy cas- 
cade also reverses but the kinetic energy added to this 
band due to wave-wave interactions is small. 

Using the same procedures used with the energy 
spectra, tests of the statistical significance of the differ- 
ences in I„ indicate greater than 99 percent confidence 
in the differences at n = 3, n = 8 and n = 13. At n 
= 3, the difference is due to a much smaller rate of 
gain of K 3 due to nonlinear interactions for LPA events. 



Fig. 3. Enstrophy E n versus n for the 500 mb level for LPA events 
(solid line); non-LPA events (dashed line) and for the 15 January 
climatological mean stationary waves (dotted line). 
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Table 2. Time mean or /„, the nonlinear kinetic energy exchange, for the various data samples. 

In this and subsequent tables, numbers to the right of the braces indicate sums over each wavenumber group. Units arc 10~ T m 1 s'*. 


All non-LPA Atlantic LPA Pacific LPA 

» days All LPA days days days 
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-180 
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-138 


-56 

13 
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-169 
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-142 
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15 
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16 
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-26 
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-72 

18 

-25 


-38 
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26 
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22 


35 


31 


35 
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26 


17 


13 


15 

30 

19 




<2 J 


36 


The difference at n - 8 represents a large increase in 
the rate of gain of A* due to nonlinear interactions and 
the difference at n = 13 represents a very large increase 
in the export of kinetic energy from n * 1 3 during LPA 
events. These statistically significant changes are con- 
sistent with the observed differences in the respective 
wavenumber ensembles in which these values of n oc- 
cur (Table 2). These changes can easily be interpreted 
to indicate that on average during LPA events a greater 
amount of K m (or £*) is exported from the transient, 
high-index scales of motion (specifically n = 13) to the 
n = 8 region of the spectrum with less cascade toward 
very large scales, especially n * 3. 

The large, significant changes in /, at n * 3 and n 
* 13 are not accompanied by significant changes in 
K n at these index values (although Kj declines as one 
might expect if I m is an important term in its budget), 
suggesting that nonlinear effects are compensated in 
some way. However, at n * 8 the significant increase 
in K„ is accompanied by a significant increase in the 
nonlinear contribution toward its maintenance. The 
increase in A, at n = 7, on the other hand, requires an 
alternative, most likely baroclinic, energy source. 

Our results by no means afford an explanation of 
the existence of the LPAs, but they do indicate the 


potential importance of cyclone-scale and small-scale 
transients in the maintenance of the LPA patterns, at 
least in a statistical sense. Saltzman (1959) first sug- 
gested that wave-wave interactions could act to main- 
tain the quasi-permanent centers of action. More re- 
cently, Holopainen and Oort (1981) revived this idea 
and showed that transient eddies act to maintain the 
local energy and enstrophy balance of midlatitude, 
time-mean, standing waves as noted earlier. In the 
present context, wave-wave interactions act to dimin- 
ish the K h and E„ for modes in the 6 < n < 9 bands 
during “normal” circulations. When the LPA subset 
(roughly '/j of the total sample) is examined, the noted 
reversal of this tendency appears. 

We can also construct results for interactions in the 
one-dimensional index m, the zonal wavenumber, but 
these are not directly comparable to conventional 
studies of cascading in this one-dimensional index be- 
cause the present results include interactions amongst 
all wavenumbers including the zonal mean wind com- 
ponents, m = 0. Conventional studies treat this wave- 
mean flow interaction separately (Saltzman, 1970). 
Nonetheless, if we compute the difference between 

A, = 2 h m or J m =X J” 

n n 


( 15 ) 
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Table 3. Time mean of J „ , the nonlinear enstrophy exchange, for the various data samples, as in Table 2. Units are 10' 1 * s'*. 
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for the LPA events and the remainder of the dataset, 
differences similar to those found by Hansen and Sutera 
(1984) appear for the Atlantic cases (Table 4). In the 



Fig. 4. The difference between the kinetic energy c a scade for LPA 
events and all other times, a/., as a function of n. Results are only 
plotted for n < 20. For higher values of n, the differences are small 
(See Table 2). See text for interpretation of this figure. 


mean for these cases, zonal wavenumbers 1-4 gain a 
large positive increment in K m and E m due to barotropic 
processes while the zonal mean flow (m = 0) receives 
less E and substantially less K during the Atlantic LPA 
events. For the Pacific cases, the zonal mean flow ex- 
periences even greater deficits from barotropic pro- 
cesses, but the enhanced upscale transfer from inter- 
mediate to large-scale waves does not appear. 

In addition, results for a meridional index, n — m, 
can be constructed. We have arbitrarily grouped these 
waves in ensembles for n - m - 1 to 5, n - m = 7 to 
15, and n - m * 17 to 29 for comparisons (Table 4 
and Table 5). Note that the major differences between 
LPAs and the mean for other days in the sample occur 
at the low wavenumber end of the spectrum at n - m 
= 5. Modes with n — m = 5 roughly double their gain 
in kinetic energy, and their enstrophy cascade reverses. 
Similar characteristics appear when the LPA events 
are subdivided into Atlantic and Pacific cases (Table 
5). Notice that the general characteristics of the non- 
linear cascading in the meridional index are very sim- 
ilar to studies of cascading in the more conventional 
zonal index, m (Steinberg et al., 1971). In general, the 
differences in the one-dimensional indices at 1 « m 
< 4 and n - m m 5 are consistent with the two-di- 
mensional changes for 6 < n < 9. The cascade of energy 







I March 1986 


BROWN, JR., HANSEN AND PANDOLFO 


483 


Table 4. Diff e rences, A lj and A Jj, in the nonlinear exchange of energy and enstrophy between the noted LPA events and the normal 
statistics for spectral index j. The results are listed for the various wave ensembles. Units for AT, are 10~ T m* s~* and for AT, ate IQ' 1 * s~ J . 



All LPA 

Atlantic 

Pacific 

All LPA 

Atlantic 

Pacific 


cases 

LPA cases 

LPA cases 

cases 

LPA cases 

LPA cases 
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A J. 
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—46 
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2 
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67 
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-223 
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-116 
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A A. 
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-74 
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-24 
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-35 
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46 

158 

-3 
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61 
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5-12 

69 
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13-30 

25 

36 

18 

40 

48 

31 

n-m 


a/ m 





1-5 

325 

237 

318 

79 

31 

91 

7-15 

-372 

-297 

-351 

-136 

-108 

-122 

17-29 

47 

60 

33 

57 

-77 

31 


into the portion of the spectrum representing LPAs 
comes from eddies of smaller zonal and meridional 
scale. 

5. Conclusion 

Locally defined height anomalies leave statistically 
significant signals in hemispheric-scale circulation sta- 
tistics. The major average differences in the kinetic en- 
ergy and enstrophy spectra between the large, positive, 
persistent S00 mb height anomalies identified by Char- 
ney et al. (1981) and more normal circulations occur 
in the two-dimensional spectral index band 6 < n 
< 9. The LPA events achieve both greater K„ and £, 
on average in this region of the spectrum. Differences 


elsewhere are negligible. Statistically significant differ- 
ences in the nonlinear cascading occur at n = 3, where 
less K and £ are gained; at n = 8, where a tripling in 
the amount of K and £ gained occurs; and at n = 13, 
where a very large increase in the amount of K and £ 
lost occurs. Only the change in cascading at n = 8 is 
accompanied by a statistically significant change in K n 
or £„. These results imply that during LPA events, 
smaller two-dimensional scale transients add K. and E 
to the wavenumbers characterizing the LPAs whereas 
during more zonal circulations wave-wave interactions 
act to dissipate this ensemble. Examination of differ- 
ences in the cascading in 1 dimensional indices indi- 
cates that the transport of energy and enstrophy to the 
LPAs comes from eddies of both smaller zonal and 


Table 5. Time mean of 7_ - 2 1" and JU ■ X A", the nonlinear exchanges of K and E in the meridional index, 

n — m. Wavenumber group sums are given to the right of braces. Units for /„ are IQ' 1 m* s~* and for ■/„ are I0~‘* s' 1 . 
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meridional extent. This suggests that the transients may 
play a regime dependent role in the dynamics of per- 
sistent circulation anomalies. The relative importance 
of these barotropic processes is an important question 
for further study. The study should be extended to de- 
termine the role of transient eddy heat transports in 
the dynamics of the anomalies. These results certainly 
do not imply that turbulent interactions alone are re- 
sponsible for the existence of the anomalies. However, 
they do imply that the transients cannot be ignored or 
treated strictly as dissipative agents. 
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APPENDIX 

Estimate of Truncation Error 

In order to examine the energetics of LPAs and of 
normal flow, it is necessary to rely upon results that 


depend on derivatives of the basic variables in spite of 
inherent problems in evaluating such derivatives. The 
presence of higher order derivatives in the Jacobian 
form of ( 10), along with the nonlinearity of the expres- 
sion, cause particular computational problems in eval- 
uating the q", r„ m (Eq. 13) that determine the wave- 
wave interaction coefficients. These derivatives may be 
calculated by finite differences or by term-by-term dif- 
ferentiation of series expansions. Finite differences in- 
troduce considerable truncation error, and so we have 
chosen the latter method of calculation. The termwise 
differentiation approach is not without problems either, 
due to the fact that the higher wavenumber, less reliable 
(i.e., more poorly resolved) components are amplified 
substantially in the differentiation process. Construc- 
tion of the expansion for f (=Vfy), for example, entails 
multiplication of ^-coefficients by the factor n(n + 1) 
by virtue of the relation 


JcosmX) —1 fcosmX) 

V 2 . . }P„"V) - — «(« + IK . . /V"(m). 

IsinmX ) fl tsin/nX J 


(A-l) 


It follows further that in the series for df/dp, the 4'- 
coefficients (A„ m , B„") are multiplied by factors that 
are 0(n 3 ). Derived expansions of this type are used to 
evaluate each of the for derivatives comprising the Ja- 
cobian. 

When the expansions for the derivatives are used to 
form the product terms found in (10), interactions oc- 
cur between waves of different index. In particular, high 
wavenumber harmonics interact to form lower-wave- 
number harmonics in the resultant expansion for J. 
Due to this situation, the calculated components of J 
are found to vary even at the lower wavenumbers as 



Fig. S. J, for the 24 LPA days in January 1963 for three different truncations. 
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the truncation limit is changed. Increasing the trun- 
cation limit of the individual series does not guarantee 
convergence of the results, however, since raising the 
limit serves to append components that are more 
poorly resolved and at the same time weighted more 
heavily. The high-wavenumber error in the differen- 
tiated series is distributed to the lower-wavenumber 
components as a result of the multiplications. Figure 
5 shows the behavior of J m versus n at 500 mb where 
the spectral coefficients have been averaged over the 
24 LPA days of January 1963. The computations have 
been performed using triangular truncation with limits 
1 8, 24 and 30. As the truncation is increased, the results 
for the 3 different truncations converge quite well for 
roughly n < 14. Comparing the truncations at n * 24 
and n = 30 alone, the results appear to converge rea- 
sonably well for roughly n < 18. 
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